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Introduction 

Sleep  restriction  and  deprivation  have  profound  negative  effects  on  cognitive  ability  and  task  performance 
(e.g.,  vigilance).  At  the  same  time,  physical  and  psychological  stressors  often  lead  to  sleep  disruptions,  which 
compromise  the  body’s  ability  to  reap  the  restorative  benefits  of  sleep.  Military  personnel  often  suffer  from 
decreased  quality  and  quantity  of  sleep,  increasing  their  susceptibility  to  a  host  of  neurological  problems  and 
limiting  their  ability  to  perform  the  challenging  mental  tasks  that  their  missions  require.  Further,  the  effects  of 
pharmacological  interventions  aimed  at  ameliorating  the  deleterious  effects  of  both  sleep  loss  and  mental  and 
physical  stress  are  only  poorly  understood  and  may  result  in  unanticipated  long-term  effects  in  those  serving 
under  combat  conditions.  Advances  in  imaging  technology  have  enabled  human  studies  of  neurochemistry, 
energy  metabolism  and  functional  brain  networks  that  were  previously  impossible. 

Our  research  uses  a  combination  of  human  imaging-based  studies  to  advance  our  fundamental  knowledge  of 
the  effects  of  sleep  and  sleep-related  stressors  on  neuroperformance.  In  order  to  achieve  this  goal,  the  team 
utilizes  ultra-high  resolution  imaging  hardware,  7T  MRI  scanner  and  High  Resolution  PET  (HRRT)  scanner  in 
tandem.  These  scanners  at  Neuroscience  Research  Institute  at  Gachon  University  are  conveniently  positioned  so 
to  be  linked  by  a  patient  tray  system,  unique  around  the  world,  that  can  carry  the  subject  from  one  scanner  to  the 
other,  thereby  keeping  the  reference  frames  of  each  modality  intact  and  physically  registered.  The  removal  of 
software  image  registration  from  data  analysis,  as  well  as  prevention  of  signal  loss  by  the  absence  of 
physical/electrical  interactions  between  the  devices  results  in  unprecedented  high  quality,  ultra-high  resolution 
images. 

Adding  this  unique  high  precision  setup  to  our  capabilities,  the  scope  of  this  research  is  to  1)  use  PET/MRI 
studies  to  understand  glucose  metabolism  and  dopamine  binding  in  the  brain,  particularly  in  the  brainstem  and 
thalamus;  2)  EEG/fMRI  studies  focusing  on  the  functional  connectivity  of  networks  between  the  thalamus  and 
cortex  that  control  the  descent  into  sleep  and  are  altered  by  sleep  deprivation;  3)  EEG/fMRI  studies  designed  to 
improve  our  understanding  of  the  interaction  between  sleep  loss,  emotional  stress  and  cognitive  function,  and 
4)  Development  of  reference  image  databases  for  high  field  MRI  studies  of  the  brain. 

Body 

The  research  accomplishments  associated  with  each  task  outlined  in  the  approved  Statement  of  Work  for  year  1 
are  as  follows: 

A.  PET/MRI  studies  of  brainstem  and  cortex  sleep  regulation 

Al.  Provide  sleep  neuroscience  expertise  to  the  Neuroscience  Research  Institute  (NRI)  of  Gachon 
University  for  PET/MRI  studies. 

Al.l.  Initial  training  in  the  Neurobiology  of  Sleep  and  Sleep  Medicine  was  provided  by  Dr.  Larson-Prior 
beginning  in  November  of  2010.  Training  included  both  didactic  information  and  hands-on  practical 
training  in  sensor  application  and  quality  control,  reading  of  the  acquired  electroencephalographic  (EEG) 
signal,  and  subject  comfort.  During  this  pre-award  training  period,  preliminary  studies  provided  important 
information  to  study  design  and  implementation. 

A1.2  Workshop  on  Study  of  Neuro-Stress  with  MRI  Tractography  and  PET  Imaging 
To  achieve  an  efficient  transfer  of  expertise  with  a  wide  foundation,  we  organized  an  International 
Workshop,  titled  Study  of  Neuro-Stress  with  MRI  Tractography  and  PET  Imaging  at  Neuroscience 
Research  Institute,  Gacheon  University,  Inchon,  Korea  Feb  16-17,  201 1.  Invited  speakers  were  specifically 
chosen  from  United  States  and  Korea  to  facilitate  the  exchange  of  ideas  between  the  researchers.  The 
workshop  included  presentations  from  different  disciplines,  PET-MR  fusion  for  neuroimaging,  animal 
models  of  traumatic  injuries,  fundamentals  of  diffusion  imaging,  application  of  multi-modality  imaging  in 
sleep  studies  and  an  introduction  to  prospects  of  connectomics  (see  Appendix  1  for  the  participant  and  talk 
lists).  Presented  below  is  the  workshop  announcement: 


Recent  advances  in  diffusion  weighted  (DW)  magnetic  resonance  imaging  (MRI)  provided  non-invasive,  in 
vivo  means  for  assessment  of  white  matter  integrity  in  human  subjects.  The  state  of  the  art  instrumentation 
based  on  7T  MRI,  high  resolution  PET  imaging  and  novel  mathematical  approaches  of  data  analysis  and 
visualization  made  it  possible  to  envision  the  use  of  DWI  to  study  white  matter  integrity  along  with  nerve 
fiber  connections. 
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Consequently,  recent  studies  of  several  diseases  and  disorders  unveiled  accurate  information  about  nerve 
fiber  pathways  and  white  matter  integrity  using  DW-MRI.  This  provided  improved  guidance  for  diagnosis, 
prognosis  and  treatment  of  acute  and  chronic  symptoms  in  neurological  diseases  and  psychiatric  disorders 
such  as  stroke,  multiple  sclerosis,  epilepsy,  Alzheimer’s  disease,  Parkinson’s  disease,  schizophrenia  and 
depression  and  traumatic  brain  injury  (TBI). 

In  addition  to  these  disorders,  the  common  and  significant  problem  of  stress  itself  has  been  associated  with  a 
number  of  physiological  and  psychological  problems  including,  but  not  limited  to,  anxiety,  depression, 
chronic  sleep  loss/insomnia,  hypertension,  migraine  and  other  forms  of  headache,  and  poor  decision  making 
(executive  functioning).  The  response  to  acute  stress  is  likely  to  affect  functional  brain  connectivity  without 
producing  significant  changes  in  anatomical  connectivity.  However,  the  effects  of  chronic  stress  (continuous 
over  more  than  1  year  as  encountered  by  warfighters)  are  unknown.  For  warfighters  both  on  active  duty  and 
home  after  having  served  under  these  conditions,  there  is  ample  evidence  of  chronic  stress  responses  that  are 
not  alleviated  simply  by  removing  the  individual  from  the  field. 

The  successful  developments  of  diffusion  weighted  neuroimaging  raise  the  possibility  of  expanding  the 
horizons  of  DW-MRI  towards  an  assessment  of  stress  related  effects  with  the  aim  of  enabling  stress 
resistance.  This  workshop  will  discuss  the  potential  of  macro-connectomics  inferred  from  tractography 
methods  as  an  aid  for  the  identification  of  those  individuals  at  greatest  risk  for  deleterious  outcomes,  for  the 
development  of  effective  treatment  strategies  and  in  the  evaluation  of  the  efficacy  of  interventions. 
Additional  discussion  will  focus  on  the  importance  of  multi-modal  neuroimaging  and  functional  imaging  in 
determining  the  interactions  between  anatomical  connectivity  and  pathologies,  including  the  role  of  high 
resolution  MRI,  functional  MRI  (fMRI),  high  resolution  molecular  imaging  using  HRRT  (High  Resolution 
Research  Tomography),  positron  emission  tomography  (PET)  and  electroencephalography  (EEG). 

This  workshop  is  planned  at  the  Neuroscience  Research  Institute  of  Gacheon  University,  Incheon,  Korea, 
home  of  state  of  art  neuroimaging  instruments  including  7T  human  MRI,  High  Resolution  PET(HRRT),  3T 
MRI  and  ZOOM  PET.  The  institute  maintains  its  cutting  edge  status  with  a  new  project  for  the  development 
of  a  14T  MRI  system. 


Workshop  Feedback:  Plenary  speaker  LG.  (ret).  James  Peake  indicated  after  praising  the  content  of  the 
workshop  that  the  knowledge  that  will  be  obtained  will  outreach  beyond  the  scope  of  sleep  studies  and  will 
help  us  better  understand  and  remedy  the  PTSD-TBI  continuum  which  will  affect  a  large  population  of 
returning  warfighters.  COL.  Karl  Friedl  emphasized  the  importance  of  this  opportunity  which  allows  the 
convergence  of  different  fields,  medicine,  neuroscience,  physics  and  engineering  in  order  to  achieve 
breakthrough  developments  that  will  address  not  only  battle  field  originating  injuries  such  as  TBI  and  PTSD 
but  also  accelerated  neurodegenerative  diseases,  e.g.  Alzheimer’s  disease,  that  is  expected  to  affect  a  large 
portion  of  the  aging  population  in  the  near  future. 

The  workshop  included  presentations  from  different  disciplines,  PET-MR  fusion  for  neuro imaging,  animal 
models  of  traumatic  injuries,  fundamentals  of  diffusion  imaging,  application  of  multi-modality  imaging  in 
sleep  studies  and  an  introduction  to  prospects  of  connectomics  (please  refer  to  the  program  in  the  Appendix 
lfor  the  participant  and  talk  lists). 
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Double  Injection  Method 


^Ippnirm 


Figure  1.  The  initial  study  design  called  for  two 
injections  under  the  assumption  that  the  radioactive 
decay  curve  would  remain  constant  between  injection 
periods.  This  allowed  subtraction  of  residual 
radioactivity  fromthe  first  injection  to  determine  the 
correct  values  for  responses  to  the  second  injection. 


A2.  Analysis  of  pre-existing,  de-identified  PET/MRI  data. 

A2.1.  We  are  currently  working  with  the  IRB  to  obtain  approval  for  use  of  existing  data  archives  at  Gachon 
University.  In  the  meantime,  we  have  modeled  and  designed  an  optimal  protocol  for  future  PET  data  studies 

that  will  be  provided  by  Gachon  team.  Our 
initial  study  design  provided  for  injection  of 
FDG  at  two  points  during  the  experimental 
period  of  overnight  sleep  (Figure  1).  The 
implementation  of  this  design  required  that  the 
second  injection  occur  at  an  asymptotic  point  on 
the  FDG  uptake  curve  while  at  the  same  time 
providing  for  only  minor  loss  of  radioactivity  in 
our  injection  sample  by  the  second  time  point. 
To  achieve  this  2-4  hours  between  injections 
(Figure  2A)  needed  to  be  allowed.  Preliminary 
analysis  of  dual  injections  suggested  that  our 
assumption  of  an  asymptotic  concentration  in 
residual  radioactivity  following  the  first  inject 
was  not  being  met  (Figure  2B).  To  further  test 
for  decay  dynamics,  we  implemented  a  dynamic 
scan  protocol,  which  verified  that  levels  were 
not  consistent  and  precluded  an  accurate 
assessment  of  metabolic  activity  the  initial 
experimental  design  (Figure  2).  We  then  devised  to  a  protocol  that  allowed  only  a  single  FDG  injection  per 
study  (Figure  3).  This  protocol  required  that  only  one  sleep  stage  (either  slow  wave  sleep  [SWS]  or  rapid 
eye  movement  sleep  [REM])  could  be  investigated  in  any  individual  subject  during  any  one  experimental 
night.  As  our  protocol  required  one  night  of  acclimatization  sleep,  these  data  were  used  to  make  a 
preliminary  decision  as  to  whether  the  Gachon  team  would  provide  data  in  REM  or  SWS  in  that  subject. 

This  was  necessary  as  we  required  30  minutes  of 
in-stage  uptake  of  our  radio-nucleotide  at  different 
times  during  a  normal  night  of  sleep  (SWS  periods 
are  longest  early  in  the  night  while  REM  periods 
are  longest  late  in  the  night).  Gachon  team 
successfully  tested  this  protocol  and  are  confident 
that  they  can  use  it  for  the  planned  studies. 

As  planning  for  the  proposed  PET  and  MRI 
studies,  Gachon  team  also  examined  typical  sleep 
patterns  in  potential  volunteers.  Questions  were 
asked  about  the  number  of  hours/night  they  slept, 
their  normal  bedtimes  and  wake  times,  and  an 
overall  score  for  sleep  quality  based  on  the  Korean 
version  of  the  Pittsburgh  Sleep  Quality  Index 
(PSQI).  The  Epworth  Sleepiness  Scale  provided 
information  on  daytime  sleepiness  in  general, 
while  the  Stanford  Sleep  Index  was  used  to  assess 
each  subject’s  level  of  sleepiness  at  a  particular. 

For  the  entire  study  group,  subjects  slept  an 
average  of  6.9±1.0  hours  per  night.  Most  subjects 
reported  no  sleep  disturbances  (PSQI  =  5.2±1.8; 
normal  sleep  reports  PSQI  <10),  although  one 
subject  reported  mild  sleep  disturbance  with  a 


Figure  2.  The  assumed  response 
required  for  accurate  calculation  of 
the  response  to  second  injection  is 
shown  in  A.  In  B  was  the  true 
response,  indicating  that  we  were 
not  recovering  an  accurate 
response  to  the  second  injection 


PSQI  of  10  (a  PSQI  of  >21  is  considered 
abnormal).  The  range  of  reported  normal  bedtimes  was  1 1PM  -  4AM,  with  most  subjects  reporting  normal 
bedtimes  around  2  AM.  Although  this  investigation  was  not  part  of  the  statement  of  work  and  therefore 
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outside  of  the  scope  of  this  project,  it  is 
described  here  because  it  informs  us  that  the 
potential  study  population  should  be 
appropriate  for  our  planned  imaging  studies. 

The  focus  of  these  studies  is  on  subcortical 
control  of  sleep  in  normal  human  subjects.  Our 
interests  lay  primarily  in  an  examination  of 
brainstem  structures,  which  are  known  to  play 
a  critical  role  in  the  regulation  of  the  ultradian 
sleep  cycle  -  the  cycle  which  oscillates 
between  REM  and  non-REM  sleep  states.  In  a 
normal  human  sleeper,  there  will  be  between  2 
and  7  such  cycles  over  a  full  night  of  sleep.  To 
date,  while  much  information  on  subcortical 
regulation  or  response  to  sleep  has  been 
gleaned  from  animal  studies,  there  are  no 
studies  examining  these  regions  in  human 
subjects  due  to  the  lack  of  spatial  resolution 
inherent  in  our  non-invasive  methods  for 
imaging  these  brain  regions.  Data  acquired  in 
cat  or  rodent  models  indicate  that  the  midline  raphe  nuclei  are  highly  active  in  wake,  show  reduced  activity 
in  NREM  sleep  states  and  are  turned  off  in  REM  sleep.  This  would  lead  us  to  expect  a  reduction  in  FDG 
label  as  we  moved  from  wake  to  SWS  to  REM. 

A3.  Collaborative  development  of  new  data  processing  and  analysis  methods. 

A3.1.  Diffusion  MRI: 

The  new  diffusion  MRI  method,  Complete  Fourier  Direct  MRI  (CFD-MRI),  developed  by  Dr.  Ozcan 
provides  an  accurate  way  of  measuring,  quantifying  and  characterizing  the  motion  of  the  molecules  within 
the  biological  tissue.  CFD-MRI  is  different  than  the  existing  methods  by  its  simple  presentation  as  a  Fourier 
transform  of  the  joint  distribution  function  of  the  number  of  spins  at  a  given  location  with  given 
displacement  integrals  (1-3)  (see  also  (4)  in  the  Appendix  2  for  a  detailed  description).  The  method  does  not 
rely  on  any  assumptions  on  the  stochastic  properties  (e.g.  Markovian  property)  or  structural  properties  (such 
as  being  symmetric)  of  the  motion.  This  generality  makes  CFD-MRI  a  truly  free  model.  Accordingly,  the 
method  provides  better  information  about  the  microstructure  of  the  tissue  under  investigation,  leading  to 
more  accurate  description  of  different  types  of  pathologies,  injuries  and  changes.  A  detailed  comparison  of 
the  existing  diffusion  MR  methods  has  been  provided  in  our  team’s  recent  publication  (5)  (see  the  attached 
manuscript  in  Appendix  3).  Mainly,  CFD-MRI  relies  on  the  Hermitian  symmetry  of  the  diffusion  weighted 
(DW)  MR  signal.  The  method’s  proof  concept  obtained  using  animal  scanners  and  relatively  simple  (but 
time  consuming)  pulse  gradient  spin  echo  (PGSE)  imaging  protocol  is  being  translated  into  7T  human 
scanners  to  use  with  the  fast  scanning  methods.  There  are  significant  challenges  that  arise  from  the  new 
method  and  7T  technology.  First,  the  use  of  fast  imaging  techniques,  namely  echo  planar  imaging  (EPI), 
creates  distortions  in  the  images.  This  effect  is  further  enhanced  specifically  at  7T  MR  due  to  the  high  field. 
The  correction  methods  introduced  earlier  by  the  Gachon  team  (6)  is  not  applicable  to  CFD-MRI  due  to  its 
specific  high  dimensional  Fourier  nature. 

Furthermore,  the  novelty  of  CFD-MRI  also  brings  demand  for  a  significantly  larger  amount  of  sampling  in 
the  Fourier  space  using  as  well  unusual  time  parameters  compared  to  conventional  diffusion  imaging 
protocols.  These  required  addressing  technical  difficulties  originating  from  vendor  specific  definition  of 
diffusion  gradient  vectors  in  pulse  sequence  programming.  After  our  teams  resolved  these  technical  issues, 
Gachon  team  collected  initial  CFD  data  from  a  fixed  human  brain  using  the  7T  scanner  under  the 
supervision  of  Dr.  Ozcan.  Considering  the  novelty  of  the  mathematical  model  and  the  7T  hardware, 
specifically  the  artifacts  created  by  the  ultra-high  magnetic  field,  the  data  are  promising.  However,  there  is  a 
lot  of  room  for  improvement  which  will  be  achieved  according  the  analysis  of  the  data  after  mathematical 
modeling  of  the  reasons  behind  the  loss  of  Hermitian  symmetry  during  data  collection.  We  are 
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Figure  3.  Final  study  design.  Only 
one  injection  per  study  per  subject 
is  acquired  to  ensure  accurate 
assessment  of  response  to  brain 
state 


hypothesizing  that  whereas  local  loss  of  symmetry  is  due  but  not  limited  to  susceptibility  artifacts  and 
magnetic  field  inhomogeneity,  radio  frequency  (RF)  pulse  inhomogeneity  might  be  causing  the  deviations. 
We  will  continue  further  our  experiments  with  the  inclusion  of  more  fundamental  steps  as  running  them  on 
liquid  samples  in  order  to  identify  the  source  of  artifacts.  In  parallel,  under  the  assumption  that  the  deviation 
is  not  preventable  with  hardware/parameter  adjustments  in  pulse  programming,  mathematical  algorithms  for 
phase  correction  towards  the  re-establishment  of  Hermitian  symmetry  are  being  developed. 

Dr.  Ozcan  has  shared  the  findings  of  CFD-MRI  at  NRI-UPenn  Joint  Mini  Symposium,  Neuroscience 
Research  Institute,  Gacheon  University  Inchon,  Korea  on  Sept  27-28,  2012  and  at  the  recent  MICCAI 2012 
Workshop  on  Computational  Diffusion  MRI,  the  15th  International  Conference  on  Medical  Image  Computing 
and  Computer  Assisted  Intervention,  Nice,  France. 

Furthermore,  with  the  aim  of  increasing  collaborative  efforts  and  gathering  information  from  different 
research  centers  around  the  world,  Dr.  Ozcan  is  hosting  a  research  topic  for  exploring  different  techniques 
including  but  not  limited  to  diffusion  MRI  to  explore  white  matter  integrity.  The  research  topic  (article 
submission  deadline  was  Oct.  15th  2012),  which  will  appear  on  the  open  access  Frontiers  in  Physiology, 
Integrative  Neuroscience,  has  been  specifically  designed  by  Dr.  Ozcan  to  better  understand  tractography’s 
role  in  white  matter  integrity  with  the  additional  aim  of  uncovering  the  shortcomings  of  existing  methods. 
The  topic  has  attracted  participation  from  different  institutions,  Harvard,  Duke,  Oregon  Health  Sciences, 
Vanderbild  Universities,  Univ.  of  South  Florida,  Univ.  of  Wisconsin,  Leiden  University,  The  Netherlands 
and  Aarhus  University,  Denmark.  The  description  of  the  topic  is  the  following: 

http://www.frontiersin.org/Integrative_Neuroscience/researchtopics/New_Models_of_Diffusion_Weight/775 
New  Models  of  Diffusion  Weighted  MRI  Signal  and  Alternative  Methods  for  Characterization  of  White 
Matter  Integrity 

Recent  advances  in  diffusion  weighted  (DW)  magnetic  resonance  imaging  (MRI)  provided  non-invasive,  in 
vivo  means  for  assessment  of  white  matter  integrity  in  human  subjects.  The  novel  mathematical  approaches 
of  data  analysis  and  visualization  made  possible  the  use  of  DWI  to  study  white  matter  integrity  along  with 
nerve  fiber  connections.  Accordingly,  recent  tractography  studies  of  several  diseases  and  disorders  unveiled 
information  about  nerve  fiber  pathways  and  white  matter  integrity.  This  has  significantly  improved  guidance 
potential  for  diagnosis,  prognosis  and  treatment  of  acute  and  chronic  symptoms  in  neurological  diseases  and 
psychiatric  disorders  such  as  stroke,  multiple  sclerosis,  epilepsy,  Alzheimer’s  disease,  Parkinson’s  disease, 
schizophrenia,  depression  and  traumatic  brain  injury  (TBI).  A  majority  of  tractography  methods  are  based 
on  diffusion  tensor  (DT)  model  which  is  limited  in  describing  complex  tissue  structure  such  as  the  nerve 
fiber  junctions.  In  addition,  the  constraints  of  higher  order  analysis  methods,  such  as  symmetry,  prevent 
innately  an  elaborate  description  of  biological  tissue  microstructure  due  to  information  loss/deterioration. 
These  raise  the  concern  about  the  sufficiency  and  the  accuracy  of  white  matter  integrity/pathology 
information  and  metrics  obtained  with  current  methods.  Consequently,  after  achieving  the  proof  of  utility, 
DW-MRI  methodologies  are  currently  open  for  improvements  and  refinements.  The  purpose  of  this  special 
issue  is  to  investigate  new  analysis  methods  and  new  mathematical  models  with  possible  introduction  of 
new  acquisition  techniques.  This  will  be  achieved  by  returning  to  the  first  principles  of  DW-MRI  signal 
formation  in  order  to  explore  new  horizons  that  remained  veiled  due  to  constraints  and  boundaries  imposed 
by  existing  methodologies.  Particular  emphasis  will  be  given  to  approaches  that  will  allow  unconstrained 
evidence  based  discoveries  using  DW-MRI  signal.  Biological,  synthetic  and  numerical  phantoms  created  for 
test  and  validation  purposes  are  highly  encouraged. 

In  parallel,  physiological  description  of  white  matter  integrity  disruption  and  its  reflection  on  DW-MRI 
signal  is  central  to  the  theme  of  the  special  issue. 

The  aim  is  geared  up  towards  the  identification  of  the  nerve  fiber  DW-MRI  signal,  separating  the  portion  of 
the  signal  from  other  elements  within  biological  tissue,  whereby  an  accurate  calculation  of  surrogate 
markers/descriptors  of  white  matter  integrity  would  be  achieved.  Accordingly,  comparison  and/or 
combination  of  alternative  MRI  methods  (e.g.  susceptibility  imaging)  with  DW-MRI  are  in  focus  of  interest 
and  are  highly  welcomed. 

The  research  topic  is  essential  in  understanding  the  changes  in  the  white  matter  and  MRI’s  limits  to  measure 
this  changes  including  the  cases  where  prolonged  sleep  deprivation  leading  to  neurostress.  By  identifying 
the  brain  regions  that  are  most  affected,  it  would  be  possible  to  assess  the  level  the  warfighters  are  affected 
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during  and  after  return  from  active  duty.  The  topic’s  goal  also  includes  the  detection  of  traumatic  brain 
injury  (TBI)  and  its  effects  on  brain  tissue. 

A4.  Recruit  new  subjects  for  studies  of  normal  sleep  using  FDG  and  raclopride. 

Recruitment  of  new  subjects  is  deferred  until  the  completion  of  IRB  approvals. 

B.  Studies  of  thalamo-cortical  network  function  during  sleep 

Bl.  Analysis  of  pre-existing  EEG/fMRI  data  on  N1  and  N2  sleep  to  determine  changes  in  thalamo¬ 
cortical  network  connectivity  in  the  descent  to  sleep. 

Bl.l  EEG/fMRI  studies  of  sleep  in  normal  human  subjects. 

Dr.  Larson-Prior’s  laboratory  developed  the  technique  of  simultaneous  electroencephalography  (EEG)  and 
functional  magnetic  resonance  imaging  (fMRI)  at  Washington  University  in  St.  Louis,  acquiring  data  from 
normal,  neuro logically  healthy,  young  adult  subjects  during  quiet  waking  rest  (2005-2007)  and  in  the 
descent  to  sleep  (2007-2008).  Data  were  acquired  from  a  total  of  67  right-handed  individuals  over  this 
period,  57  of  whom  participated  in  studies  designed  to  examine  brain  network  activity  during  resting  30 
subjects  provided  full  and  good  quality  data  sets  for  analysis  (7).  Based  upon  analysis  of  EEG  data,  subjects 
were  divided  into  three  groups,  1)  those  who  remained  in  an  alert  waking  state  during  scanning  (n=9),  2) 
those  who  were  drowsy,  failing  into  the  transition  state  between  wake  and  sleep  (stage  1;  n=16)  and  3)  those 
who  were  overtly  asleep  (stage  2,  n=6). 

Analysis  of  these  data  utilized  the  methods  of  blood  oxygen  level  dependent  (BOLD)  functional 
connectivity  to  define  a  set  of  15 1  functional  brain  regions  of  interest  (7)  which  were  then  evaluated  for 
each  condition  for  correlated  activity  during  the  scan  session.  The  results  of  these  studies  clearly  indicate 
changes  in  resting  state  brain  network  connectivity  in  the  transition  from  alert  wakefulness  to  sleep.  As 
anticipated,  one  change  lay  in  the  relationship 
between  the  default  mode  network  (5),  which  is 
most  active  in  non-task  related  rest  signaling 
disengagement  from  the  external  world  (P),  and 
two  networks  involved  in  attending  to  and 
interacting  with  external  stimuli  (dorsal  attention 
network,  ATTN,  (10);  Executive  control  network 
EXEC,  (11)).  During  quiet  rest,  the  DMN  is 
most  active  while  the  networks  generally  most 
active  in  response  to  environmental  stimuli 
exhibit  reduced  activity  (12);  thus,  these 
networks  are  anti-correlated  in  functional 
connectivity  analyses.  As  individuals  fall  asleep, 
the  strength  of  these  anti-correlations  was  seen  to 
be  reduced  (Figure  4)  indicating  a  reduction  in 
the  strength  of  inter-network  connectivity  in  the 
descent  to  sleep. 

During  quiet  rest,  there  are  strong 
correlations  between  mediodorsal  thalamic 
regions  and  elements  of  the  DMN  (13).  If  these 
relationships  are  removed  to  enable  an  investigation  of  other  brain  regions  associated  with  thalamic  activity, 
one  finds  correlations  to  anterior  cingulated,  supplemental  motor  area,  and  the  basal  ganglia.  Further, 
differences  in  mediodorsal  thalamic  connectivity  differ  between  eyes  closed  and  eyes  open  fixated 
conditions,  with  the  fixation  condition  exhibiting  positive  correlations  to  the  attentional  networks  (ATTN). 
During  the  descent  to  sleep,  thalamic  connectivity  changes  were  primarily  noted  within  the  thalamus,  with 
correlations  strengthening  between  5  thalamic  regions  from  wake  to  stage  2  sleep  (7). 

B2.  Data  collection  and  analysis  to  (1)  increase  the  number  of  subjects  reaching  N2  sleep,  and  (2) 
collect  data  from  subjects  reaching  N3  sleep. 

We  will  start  data  collection  upon  the  completion  of  IRB  approvals. 


■  DMN:ATTN 
□  DMN:EXEC 


-0.30  Wake  Drowsy  Sleep 

Figure  4.  Functional  connectivity  analysis  of 
subjects  in  the  descent  to  sleep  indicated  that  the 
anti-correlated  activity  between  the  default  mode 
network  (DMN)  and  two  networks  known  to  play  a 
role  in  attentional  regulation  and  cognitive  control 
(ATTN  and  EXEC)  is  significantly  reduced  from  alert 
wakefulness  to  early  sleep  stages 
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B3.  Collect  normative  data  on  cognitive  vigilance  (PVT)  prior  to  and  following  sleep  (task  flanks 
sleep)  to  analyze  thalamo-cortical  network  interactions. 

Normative  data  collection  will  be  initiated  upon  the  completion  of  IRB  approvals. 

C.  Assemble  and  annotate  a  reference  dataset  of  7T  brain  images. 

The  7T  reference  data  set  will  be  acquired  after  the  completion  of  IRB  approvals. 

Preparations  in  Anticipation  of  Future  Data  Collection  and  Analysis 

In  order  to  address  the  variety  of  neuroscientific  challenges  that  will  be  encountered  with  the  data  collection, 

Dr.  Kathleen  Gates,  a  neuroscientist,  joined  the  Arlington  Innovation  Center  in  August  2012.  Her  work  focuses 
on  developing,  testing,  and  applying  novel  methods  that  enable  reliable  inferences  when  studying  individuals. 

In  this  vein,  she  has  collaborated  in  research  on  mother-infant  dyads,  heart  rate  variability,  and  daily  diary 
studies  that  attend  to  person-specific  substantive  questions.  By  the  incorporation  of  Dr.  Gates’  expertise  on 
developing  statistical  methods  for  analyzing  brain  imaging  (7),  the  team’s  perspective  and  horizons  has  been 
prepared  to  tackle  the  complex  challenges  of  neuroscience  of  sleep. 

The  IT  infrastructure  that  will  provide  data  storage,  analysis  using  computational  methods  and  facilitate 
communication  between  different  centers  of  the  project  in  a  secure  manner  has  been  implemented  at  the 
Arlington  Innovation  Center.  The  network  consists  of  a  firewall  and  VPN  device.  Whereas  publicly  available 
information  are  presented  at  AIC  Webserver  (Dell  PowerEdge  R310,  4  core  Intel®  Xeon®  X3470,  16Gb 
memory,  4Tb  raid  5  disk  capacity),  the  collaborators  in  the  project  have  the  ability  to  access  their  data  over  a 
secure  VPN  connection  on  the  computational  server  (Dell  PowerEdge  T710,  8  core  Intel®  Xeon®  X5647, 

96Gb  memory,  14Tb  raid  5  disk  space)  that  is  protected  behind  AIC  firewall.  The  firewall  and  VPN  services  are 
provided  by  a  Cisco  ASA5520  appliance  and  Cisco  2960  switch.  A  beta  version  of  in-house  developed 
collaborative  platform  that  allows  the  collaborators  to  drop  and  annotate  data,  manuscripts,  presentations  and 
documents  have  also  been  installed  on  the  storage  server.  The  platform  will  speed  up  data  analysis  and 
development  of  new  ideas. 

Key  Research  Accomplishments 

•  Provided  initial  training  in  the  practical  and  theoretical  basis  of  electroencephalography  (EEG)  for  the 
study  of  human  sleep  to  student  researchers  at  the  NRI,  Gachon  University 

•  In  collaboration  with  colleagues  at  NRI,  established  a  mechanism  by  which  continued  training  in  scoring 
sleep  based  on  overnight  EEG  studies  could  be  maintained 

•  All  research  team  members  worked  to  implement  an  optimal  study  design  for  initial  investigations  into 
brainstem  regulation  of  human  sleep  using  FDG-PET  imaging 

•  IRB  approval  from  Washington  University  in  St.  Louis  was  obtained  for  the  prospective  studies  on 
cognition  in  sleep  on  9/21/2012 

•  Programming  of  behavioral  tasks  is  currently  being  done  at  Washington  University  in  St.  Louis 
preparatory  to  initiation  of  EEG/fMRI  studies 

•  Further  development  of  the  new  and  more  general  method  for  diffusion  weighted  magnetic  resonance 
imaging,  Complete  Fourier  Direct  MRI. 


Reportable  Outcomes 

Manuscripts,  abstracts,  presentations: 

•  Dr.  Larson-Prior  provided  didactic  training  in  the  Neurobiology  of  Sleep  to  students  and  faculty  at  the 
Neuroscience  Research  Institute  at  Gachon  University  in  July  2011  with  retroactive  support  from  this  grant. 
Presentation  entitled:  “Functional  neural  network  modulation  in  wake,  sleep  and  learning” 


•  In  July  2011,  Dr.  Larson-Prior  was  asked  to  speak  to  MCB  Korea  on  the  importance  of  sleep  in  humans 
with  retroactive  support  from  this  grant.  This  discussion  was  provided  for  the  lay  public  and  included  a 
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short  text  that  could  be  disseminated  by  MBC  to  the  public.  The  presentation  is  included  in  the 
Appendix  4. 

•  Dr.  Larson-Prior  presented  a  summary  of  research  progress  and  the  rationale  for  the  studies  supported 
under  this  contract  to  Michael  P.  Brazaitis  MD,  MA,  Program  officer  at  TATRC  entitled:  “Sleep  and 
Neuroperformance”  in  October  201 1.  At  the  same  meeting  Dr.  Ozcan  lead  a  discussion  on  the  diffusion 
weighted  imaging  magnetic  resonance  imaging’s  capabilities  to  detect  white  matter  integrity  using 
tractography  in  order  to  obtain  a  reference  set  for  healthy  population. 

•  Dr.  Larson-Prior  and  Dr.  Son  presented  their  work,  entitled  “Brainstem  Regulation  of  Sleep  and 
Waking:  A  PET/MRI  Perspective”  at  the  Study  of  Neuro-Stress  with  MRI  Tractography  and  PET 
Imaging  Workshop,  Neuroscience  Research  Institute,  Gacheon  University  Inchon,  Korea,  Feb  16-17, 
2012 

•  Larson-Prior,  LJ,  Power,  JD,  Vincent  JL,  Nlan,  TS,  Coalson,  RS,  Zempel,  J,  Snyder  AZ,  Schlaggar  BL, 
Raichle,  ME  and  Petersen  SE  (201 1)  Modulation  of  the  brain’s  functional  network  architecture  in  the 
transition  from  wake  to  sleep.  Prog  Br  Res,  193:277-294.  PMID  21854969 

•  McAvoy,  M,  Larson-Prior,  L,  Ludwikow,  M,  Zhang,  D,  Snyder  AZ,  Gusnard  DL,  Raichle  ME,  d- 
Avossa  G  (2012)  Dissociated  mean  and  functional  connectivity  BOLD  signals  in  visual  cortex  during 
eves  closed  and  fixation.  J.  Neurophysiol.  In  press.  PMID  22875902 

•  Gates,  K.  M.  Identifying  subgroups  using  fMRT  connectivity  maps.  Annual  Conference  for  the  Society 
for  Neuroscience,  New  Orleans,  USA,  October  13-18,  2012 

•  Molenaar,  P.  C.  M.  and  Gates,  K.  M.  A  high  precision  method  to  estimate  effective  connectivity 
networks  at  the  group  and  individual  levels.  World  Congress  for  the  International  Organization  of 
Psychophysiology,  Pisa,  Italy,  September  13-17,  2012 

•  Alpay  Ozcan,  Complete  Fourier  Direct  Magnetic  Resonance  Imaging  (CFD-MRl)  Equations  for  Diffusion 
MRI  Signal,  CDMRT12,  MICCAI  2012  Workshop  on  Computational  Diffusion  MRI,  the  15th  International 
Conference  on  Medical  Image  Computing  and  Computer  Assisted  Intervention,  Nice,  France,  1-5  October 
2012. 

•  Alpay  Ozcan,  Kenneth  H.  Wong,  Linda  Larson-Prior,  Zang-Hee  Cho,  Seong  K.  Mun:  Background  and 
Mathematical  Analysis  of  Diffusion  MRI  Methods.  International  Journal  of  Imaging  Systems  and 
Technology,  Feb.  2012,  Vol.  22,  pp.  44-52 

•  Alpay  Ozcan,  James  D.  Quirk,  Yong  Wang,  Qing  Wang,  Peng  Sun,  William  M.  Spees  and  Sheng— Kwei 
Song,  The  Validation  of  Complete  Fourier  Direct  MR  Method  for  Diffusion  MRI  via  Biological  and 

Numerical  Phantoms,  pp. 3756-3759,  33rd  Annual  International  Conference  of  the  IEEE  EMBS,  Boston, 
Massachusetts  USA,  August  30  -  September  3,  201 1 

•  Dr.  Ozcan  has  initiated  the  research  topic  in  Frontiers  in  Physiology,  Integrative  Neuroscience:  New 
Models  of  Diffusion  Weighted  MRI  Signal  for  Characterization  of  White  Matter  Integrity 

•  Dr.  Ozcan  was  the  invited  speaker  at  McLean  Hospital,  Harvard  Medical  School,  Belmont  MA,  Neuroscience 
Seminar  Series  in  June  2012.  He  presented  white  matter  integrity  and  its  measurement  with  the  Complete 
Fourier  Direct  MRI. 

•  Dr.  Ozcan  was  an  invited  speaker  at  Instituto  San  Lazaro  de  Neurociencias’  Inaugural  Symposium,  Circulo 
Medico  Salta,  Argentina,  March  28-29,  2012.  His  presentation  was  on  diffusion  MR’s  limitations  in 
measuring  the  white  matter  integrity  with  specifics  to  neurological  disorders. 

•  Dr.  Ozcan  was  an  invited  speaker  at  NRI-UPenn  Joint  Mini  Symposium,  Neuroscience  Research  Institute, 
Gacheon  University  Inchon,  Korea  on  Sept  27-28,  2011  where  he  presented  the  new  diffusion  MRI  model: 
Complete  Fourier  Direct  MRI. 

•  Workshop:  Study  of  Neuro-Stress  with  MRI  Tractography  and  PET  Imaging,  Neuroscience  Research 
Institute,  Gacheon  University  Inchon,  Korea,  Feb  16-17,  2012 
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Conclusion 

In  anticipation  of  upcoming  data  from  Gachon  team,  we  have  designed  an  appropriate  PET  imaging  protocol 
that  will  allow  us  to  differentiate  regional  metabolic  activity  between  different  sleep  cycle  periods.  Clarification 
of  timing  methodology  for  double  injection  PET  studies  constitutes  an  important  step  towards  accurate  data 
acquisition,  analysis  and  interpretation  of  results.  As  the  time  period  between  sleep  cycles  creates  inter-  and 
intra-  subject  variability,  the  variable  decay  amount  of  the  contrast  agent  can  potentially  lead  to  inconclusive 
data  acquisition.  The  establishment  of  the  protocol  will  lead  to  more  accurate  and  precise  analysis  and 
interpretation  of  future  findings  after  the  start  of  data  collection  awaiting  IRB  approval  completion. 

Pilot  diffusion  MRI  data  have  been  collected  using  a  fixed  human  brain  in  order  to  apply  the  new  and  more 
accurate  mathematical  model,  CFD-MRI,  to  MR  acquisition  for  tractography.  The  analysis  indicates  the 
necessity  of  further  investigations  for  better  understanding  and  quantifying  various  artifacts  in  data  collection. 
The  information  and  deeper  understanding  of  anatomy  and  white  matter  connectivity  will  lay  the  foundation  of 
metabolic  and  dynamic  functional  changes  within  the  brain.  The  microstructural  form  described  by  the  diffusion 
MRI  is  necessary  for  a  complete  understanding  of  brain  function.  Therefore,  the  improvements  towards  a  higher 
accuracy  of  diffusion  MRI,  rely  on  understanding  and  eliminating  the  artifacts,  perturbation  errors  and  their 
propagation  in  diffusion  MR  acquisition.  CFD-MRI  is  tackling  these  important  issues  by  going  back  to  the 
origins  of  signal  formation  and  addressing  the  fundamental  modeling  challenges. 

Using  the  simultaneous  fMRI-EEG  techniques,  brain  regions  of  interest  and  their  activity  has  been  identified 
The  results  of  these  studies  clearly  indicate  changes  in  resting  state  brain  network  connectivity  in  the  transition 
from  alert  wakefulness  to  sleep.  An  important  change  lay  in  the  relationship  between  the  default  mode  network 
and  two  networks  involved  in  attending  to  and  interacting  with  external  stimuli.  Further  data  collection  will 
provide  us  with  better  description  of  the  relationship  between  the  default  network  and  the  peripheral  networks 
that  are  activated  during  sleep  cycles. 


“So  What”  Section: 

As  the  plenary  speakers  of  the  workshop  have  indicated,  the  outcomes  of  this  Neuroperformance  project,  which 
is  built  on  convergence  of  different  disciplines,  will  have  significant  implications  not  only  for  warfighters  but 
also  the  general  population.  Specifically  our  understanding  of  form  via  high  resolution  anatomical  images  and 
diffusion  MRI  tractography  (with  information  on  white  matter  integrity)  using  7T  MRI  scanners  and  our 
understanding  of  function,  including  metabolic  activity  (via  PET),  BOLD  activity  (via  fMRI)  and  electrical 
activity  (via  EEG),  are  essential  to  a  neurobiologically  informed  approach  to  improvements  in 
neuroperformance  under  conditions  of  restricted  sleep  and  high  stress.  In  addition  to  improving  our 
understanding  of  the  effects  of  stress  and  sleep  loss  on  performance,  and  to  providing  a  firmer  understanding  of 
the  effects  of  pharmacological  interventions  under  these  conditions,  these  studies  will  provide  novel  data  on  the 
multifactorial  impact  of  sleep  loss,  stress  and  cognitive  load  in  the  development  of  neuropsychiatric  disorders 
commonly  associated  with  abnormalities  of  sleep  such  as  depression  and  anxiety. 
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Abstract.  By  carefully  re-examining  the  first  principles  of  diffusion 
weighted  magnetic  resonance  imaging  (DW-MRI)  signal  formation  and 
deriving  the  mathematical  model  from  scratch,  the  foundation  for  an  ac¬ 
curate  and  unifying  Fourier  based  theory  is  constructed.  The  derivations 
are  specifically  obtained  for  DW-MRI  signal  by  including  imaging  gradi¬ 
ents  and  without  taking  the  magnitude  of  the  signal  using  particle  meth¬ 
ods  instead  of  the  conventional  partial  differential  equations  approach. 
The  signal  is  shown  to  be  the  Fourier  transform  of  the  joint  distribution 
of  number  of  the  magnetic  moments  (at  a  given  location  at  the  initial 
time)  and  their  displacement  integrals  which  augments  the  k- space  by 
three  more  dimensions.  The  joint  distribution  function  is  recovered  by 
applying  the  Fourier  transform  to  the  complete  high-dimensional  data 
set.  A  physically  meaningful  real  valued  distribution  function  is  obtained 
after  phase  corrections  are  applied  for  the  re-establishment  of  Hermitian 
symmetry  in  the  signal.  Consequently,  the  method  is  fully  unconstrained 
and  directly  presents  the  distribution  of  displacement  integrals  without 
any  assumptions  such  as  symmetry  or  Markovian  property.  The  joint 
distribution  function  is  visualized  with  isosurfaces  overlaid  on  the  distri¬ 
bution  map  of  the  number  of  magnetic  moments  with  low  mobility.  The 
presentation  provides  a  highly  detailed  level  of  microstructural  informa¬ 
tion.  The  results  are  demonstrated  on  the  experimental  data  obtained 
from  an  ex-vivo  baboon  brain. 

Keywords:  magnetic  resonance  imaging,  diffusion  weighted  imaging, 
Fourier  transform 


1  Introduction 

Diffusion  weighted  (DW)  MRI’s  ability  to  estimate  molecular  displacements  pro¬ 
vides  information  about  the  microstructure  within  which  the  motion  takes  place. 
The  technique  has  successfully  been  used  in  a  variety  of  clinical  and  research 
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areas  such  as  cancer  diagnosis  [1],  white  matter  integrity  assessment  [2],  in  mon¬ 
itoring  of  neurological  diseases  such  as  multiple  sclerosis  [3]  and  disorders  [4]  as 
well  as  neonatal  development  [5]. 

However,  the  existence  of  several  methods  for  modeling  the  DW-MRI  signal 
is  proof  that  the  field  has  not  yet  reached  a  satisfactory  consensus.  A  thorough 
mathematical  analysis  and  a  description  of  the  subject  and  the  background  of 
existing  methods  is  presented  in  [6].  Herein,  to  provide  a  unified  approach,  the 
MRI  signal  in  the  presence  of  diffusion  sensitizing  gradients  is  derived  by  care¬ 
fully  re-examining  the  first  principles  of  DW-MRI  signal  formation  and  deriving 
the  mathematical  model  from  scratch  using  particle  methods  in  the  spirit  of  the 
work  of  McCall  et  al.  [7]  rather  than  starting  from  the  conventional  partial  dif¬ 
ferential  equations  (PDE)  of  Stejskal  and  Tanner  [8].  The  outcome  is  the  Com¬ 
plete  Fourier  Direct  (CFD)  MRI  method  [9]  that  establishes  a  high  dimensional 
Fourier  relationship  between  the  signal  and  the  joint  distribution  of  position  and 
displacement  integrals  of  the  particles. 

2  Signal  Formation 

The  derivation  of  the  CFD-MR  signal  model  is  accomplished  by  simply  express¬ 
ing  the  magnetization  of  each  magnetic  moment.  For  this  purpose,  the  effect 
of  spin-spin  relaxation  can  be  neglected.  The  evolution  of  the  transverse  mag¬ 
netization  of  the  ith  magnetic  moment  is  described  in  a  standard  manner  by  a 
rotating  magnetization  vector: 


rriiit)  =  exp(-j7i7i)mi(io).  (1) 

Here,  7  is  the  gyromagnetic  ratio,  the  transverse  magnetization  vector,  m^,  is 
written  in  complex  number  form  with  777  (to)  denoting  the  initial  magnetization 
tipped  to  the  transverse  plane, 


J2i  = 


•  Xi(r  )dr 


(2) 


describes  the  phase  as  a  function  of  the  magnetic  field  gradients  G(x,t)  G  M 3, 
and  the  position  of  the  magnetic  moment  is  27  G  JR3. 

The  signal  originates  from  all  of  the  magnetic  moments 


M(t)  =  7^  rrij  (t)  =  7exP 

i  i 


(3) 


In  order  to  precisely  obtain  the  CFD-MR  signal,  the  time  evolution  of  the  phase, 
Qi,  in  Eq.  2  must  be  explicitly  written  on  the  right  hand  side  of  Eq.  3  by 
expressing  27.  The  time  dependent  position  of  each  magnetic  moment  is  given 
as: 

Xi(t)  =  Xi(t0)  +  Wi{t)  (4) 
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where  Wi(t)  G  JR3  represents  the  displacement  of  the  magnetic  moment  from 
its  initial  position  with  Wi(to)  =  0.  The  function  Wi(t)  could  express  any  kind 
of  displacement  such  as  Brownian  motion,  molecular  movement  in  biological 
tissue  with  different  medium  and  obstacles,  coherent  motion  or  any  combination 
thereof.  The  function  does  not  have  to  be  differentiable  but  it  must  be  continuous 
to  make  physical  sense  since  a  magnetic  moment  cannot  disappear  at  a  given 
point  and  reappear  at  another.  No  other  requirements  and/or  assumptions  about 
motion  properties  in  the  time  domain,  deterministic  and/or  stochastic,  such  as 
being  Markovian,  is  required  in  the  derivations  of  the  CFD  model.  Based  on 
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Fig.  1.  The  pulsed-gradient  spin-echo  (PGSE)  pulse  sequence  and  the  definition  of  the 
variables  used  in  the  calculations.  RO  is  for  read  out,  PE  for  phase  encode,  SS  for  is 
the  slice  select  gradient,  Diff  marks  the  motion  sensitizing  pulses  and  ACQ  stands  for 
digital  acquisition  period.  In  practice,  the  MR  pulse  sequences  implement  the  rewind 
(rw)  gradients  such  that  the  amplifiers  are  turned  on  and  off  at  the  same  times.  The 
magnetic  field  gradient  pulses  during  the  radio  frequency  (RF)  pulses  are  not  shown 
because  their  effect  is  overwhelmed  by  the  RF  pulse  (see  also  the  text). 


this  model,  the  phase  of  the  ith  magnetic  moment  during  the  digital  acquisition 
period  of  the  pulsed-gradient  spin-echo  (PGSE)  sequence  of  Fig.  1  (for  /mr  =  2 
dimensional  imaging  slice)  is  obtained  after  tedious  but  routine  derivations  as 
follows.  In  the  beginning  of  each  acquisition,  a  tt/2  radio  frequency  (RF)  pulse 
tips  the  longitudinal  magnetization  to  the  transversal  plane.  The  initial  time,  to, 
is  chosen  accordingly  to  be  the  end  time  of  the  RF  pulse  when  the  magnetization 
is  fully  tipped  to  the  transversal  plane.  Using  the  definitions  of  the  variables  in 
Fig.  1,  the  evolution  of  the  phase  is  described  as: 

1)  First,  imaging  gradients  for  read  out  rewinding,  Gro  G  JR3,  phase  encoding, 
Gpe  G  JR3  and  slice  select  Gss  G  M3  on  the  time  interval  [to,trw]  are  turned  on. 
After  these  gradients  are  applied  the  phase  and  the  magnetization  become 


hw 

f^rw  =  J (Gto(t)  +Gpe(r)  +  GSs(t))  •  Xi{r,  to)  dr 

to 

(tw  > 

Atrw  Xi(to)  +  I  Wi{r)  dr 

to  ) 

=  Z\trw  (Gro  H-  Gpe  H-  Gss)  •  2^(to) 
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(5) 


trw 

+(Gro  +  Gpe  +  Gss)  •  J  Wi(r)  dr. 

to 


Note  that  the  calculations  are  carried  out  assuming  ideal  gradient  amplifiers  re¬ 
sulting  in  rectangular  shaped  gradient  pulses  in  the  formulation.  Routine  deriva¬ 
tions  for  trapezoidal  shapes  are  omitted  for  clarity. 

2)  The  7 r  RF  pulse  between  the  diffusion  gradient  pulses,  Gd  E  I?3,  and  the 
accompanying  slice  select  gradient  provide  theoretical  sign  reversal  of  the  phase. 
The  slice  select  gradient  during  this  period  does  not  encode  magnetic  moment 
motion  into  the  signal  (thus  it  is  not  shown  in  Fig.  1).  This  is  a  direct  consequence 
of  the  RF  pulse  imposes  a  strong  magnetic  field  that  changes  all  of  the  entries  of 
the  main  magnetic  field  Bq.  In  comparison,  the  slice  select  gradient’s  modification 
of  the  third  component  of  the  Bq  vector  is  negligible.  Therefore  the  strong  (RF) 
magnetic  field  rotates  the  magnetization  vector  nq  of  each  magnetic  moment 
without  being  affected  by  diffusion. 

Since  Xi(tdk )  =  Xi(t0)  +  Wi(tdk ),  k  =  1, , . . , 4;  at  t  =  tdA 


= 


tdA  td2 

=  J  Gd  •  Xi(T,t0)d,T  -  J  Gr>  ■  Xi(r,to)dT 

td3  tdl 


=  Gd 


{tdA  -  td3 )  Xi(t0)  +  /  Wi(r)  dr 

td3  ) 

f  td2 

(td2  ~  tdi)  Xi(t0)  +  J  Wi(r)  dr 

\  tdl 

(  tdA  td2  \ 

J  Wi(r)dr-  J  Wi(r)  dr 

^td3  tdl  / 

+((^4  -  tds)  -  (td 2  -  td i))  Gd  •  Xi(t0). 


-Gd 


=  Gd 


(6) 


The  sign  change  on  i?rw  is  the  effect  of  the  ir  pulse,  which  is  incorporated  in 
the  expression  for  i?D-  If  the  diffusion  gradient  times  are  equal,  i.e.,  td 4  —  td 3  = 
td 2  —  tdi  =  S,  then  the  last  term  in  Eq.  7  is  equal  to  zero,  erasing  the  influence 
of  initial  position  from  motion  encoding  part  of  the  signal.  As  a  result,  the 
formulation  with  displacement  integrals,  f  Wi(r)  dr ,  is  a  more  meaningful  choice 
than  the  center  of  mass  (COM)  of  random  walk,  f  x(r)  dr ,  introduced  in  [10]. 
3)  The  last  part  of  the  sequence  is  where  the  data  are  collected: 


,M=/ 


Gro-{Xi(t0)+Wi(T))  dr  =  (t-tg, cq)  G 


t 

i(to)+  J 


Gro-Wi(r)dr ,  (8) 


fti(t,  Gpe,  GS| 


u  Gd)  = 
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if  ^acq  ^^rw)  Gro  '  %i(f o)  Atrw  G pe  •  Xiiff)  Atrw  Gss  •  X-iftf) 
+Gd  •  Wd  +  Gro  •  Wiacq(t)  —  (Gro  +  Gpe  +  Gss)  •  W™ 


(9) 


Equation  9  incorporates  three  integrals  of  the  displacement  Wi(t)  :  {Wf,  Wacq,  W/w) 
corresponding  to  the  displacement  integrals  for  diffusion  (d),  analog  to  digital 
conversion  acquisition  (acq)  and  initial  rewind  (rw)  gradient  time  periods  re¬ 
spectively 


tdA  td2  t  ^rw 

Wi(r)  dr ,  Wr  =  j  Wi(r)  dr. 

td3  tdl  ^acq  to 

(10) 

Moreover,  in  Eq.  9,  the  first  two  terms  that  multiply  the  initial  position  Xi(to) 
are  the  foundation  of  MR  imaging  and  can  be  written  in  the  fc-space  formalism 
as 


W?  =  j  Wi(r)  dr  -  J  Wi(r)  dr,  W^cq(t)  =  J 


kmr (t,  Gpe)  ( t  taCq  Z\trw)  Gro  Atrw  G 


pe 


(11) 


where  Gro  and  Atrw  are  fixed,  and  t  and  Gpe  define  the  sampling  points  in  k- 
space.  Also,  the  term  cps iice  =  —  Atrw  Gss  •  Xiftf)  is  constant  for  all  i  because  the 
slice  select  axis  component  of  Xi(to)  is  the  slice  position. 

Without  loss  of  generality,  the  coordinate  frame  used  for  the  equations  is 
the  imaging  coordinate  frame  defined  by  the  read-out,  phase  encode  and  slice 
select  gradients,  Gro  =  \gTOi,  0,  0],  Gpe  =  [0,  gpe 2,  0],  Gss  =  [0,  0,  gss3]-  Using 
this  notation,  Eq.  11  implies  that  t  =  kmri/groi  +  £acq  +  Atrw  and  gpe 2  = 
— kmr2/Atrw  where  kmr  has  its  components  equal  to 


[^mrl?  kmr2i  ^mr3]  [(f  ^acq  AtTW)  gYO\,  Atrw  gpe2i  0]  • 

With  this,  Wfcq(t)  becomes  a  function  of  kmri  and  it  is  denoted  by  Wfcq(kmri). 
The  coefficients  of  Wrw  in  Eq.  9  are  written  as  a  vector: 


krw  —  [  Qroli  kmr2/Atrw,  £/ss3]  •>  (12) 

which  is  an  affine  function  of  kmr. 

By  defining  kr>  =  Gd,  Eq.  9  is  rewritten  concisely: 

fti  =  kmr  •  Xi(to)  +  &d  •  +  krw  •  Wfw  +  W^q(kmr i)  gYOi  +  ^siice?  (13) 

reflecting  the  partition  of  the  effect  of  initial  position  and  displacement  integrals 
on  the  phase  of  each  magnetic  moment.  Since  (^siice  is  constant  for  all  i,  it  is  taken 
out  of  Eq.  13  with  the  appropriate  rotation  of  the  magnetization  coordinate 
frame  on  a  slice  by  slice  basis. 

In  Eq.  13,  the  pairs  (fcmr,  #i(to))>  (&d,  Wf)  and  (fcrw,  Wfw)  are  clearly  Fourier 
dual  variables  while  Wacq  introduces  dependence  on  magnetic  moment  motions. 
The  derivations  of  the  total  signal  expression  in  the  next  section  clarifies  this 
Fourier  relationship. 
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3  Complete  Fourier  Direct  MRI 


The  total  CFD-MRI  signal  comes  from  all  of  the  magnetic  moments  and  it  is 
the  sum  of  individual  magnetizations  with  the  real  valued  initial  magnetization 
rrii(to)  =  mo.  Equation  3  can  be  re-written  using  the  relationship  in  Eq.  13  using 
the  Fourier  variables  of  Section  2: 

^cfd(^mr5  ^rw)  —  ^0  E  exp(— &D,  fcrw)).  (14) 

i 

The  dimension  of  the  signal  domain  is  (/mr  +  3  +  3)  but,  whereas  (kmr,kr>)  £ 
jj>+r+3  are  independent  variables,  krw  G  JR 3  is  affinely  dependent  on  kmr  by 
Eq.  12.  Accordingly,  the  CFD  k- space  variable  is  defined  by  joining  the  inde¬ 
pendent  variables:  kc fd  =  (Aw^d)  G  1 Rl™+3. 

Since  the  properties  of  interest,  e.g.,  tissue  structure,  are  location  dependent, 
instead  of  adding  up  the  magnetization  magnetic  moment  by  magnetic  moment, 
the  sum  in  Eq.  14  can  be  carried  out  in  the  joint  position-displacement  integral 
space.  For  this  purpose,  the  total  CFD  distribution  function,  Pc^tal(x,  VF),  is  de¬ 
fined  here  as  the  number  of  magnetic  moments  with  the  initial  position  x  G  JR 3  at 
time  to,  possessing  the  displacement  integral  values  of  W  =  (VFd,  VFrw,  VFfcq)  G 
JR7 .  By  definition,  P^tal  is  real  valued  and  positive  semi-definite. 

All  of  the  individual  magnetic  moment  signals  coming  from  the  infinitesimal 
volume  element  dx  dW  around  the  (position,  displacement  integral)  point  (x,  W ) 
must  be  integrated  to  form  the  DW-MRI  signal.  The  sum  in  Eq.  14  becomes 
(absorbing  mo  into  P^tal  for  ease  of  notation): 

Scfd(fcmr,fcD,fcrw)  =  /  ^cfd'+WOx 

exp(— jy  (fcmr  •  x  +  feD  *  Wd  +  kYW  •  Wrw  +  Wfcq(fcmri)  gv Qi))  dx  dW  (15) 

defining  the  Fourier  transform  of  the  total  CFD  distribution  function  P^al 
evaluated  at  the  frequencies  specified  by  (fcmr,  /+>,  tw)  : 

S'cfd^mr,  &D,  tw)  =  mr,  &D ,  kTW)  (16) 

with  nonlinearities  added  by  lFacq.  Accordingly,  the  only  physical  constraint 
that  guides  the  signal  reconstruction  is  the  following. 

Since  P^al  is  real  valued,  Scf a  is  Hermitian  symmetric. 

Experimentally,  P^tal  cannot  be  obtained  by  the  inversion  of  Eq.  16  via  discrete 
Fourier  transform  (DFT)  because  krw  cannot  be  independently  sampled  at  nec¬ 
essary  points  (e.g.,  on  a  Cartesian  grid)  due  to  the  affine  dependence  of  kYW  on 
kmr  in  Eq.  12.  However,  P*fatal  is  not  the  main  point  of  interest  for  discovering 
the  microstructural  properties  of  the  environment.  Among  the  elements  of  W, 
the  focus  is  on  the  most  descriptive  MRI  observable,  Wd.  The  objective  is  to 
marginalize  Wrw  from  Eq.  15  to  obtain  the  distribution 

Pcm(x,Wd)  =  j  P^\x,Wd,  W™)dW™=  7+rjfcD){5cfd(fcmr,fcD,0)}.  (17) 
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without  being  able  to  sample  at  kYW  =  0. 

With  this  issue  at  hand  and  the  disruption  of  the  Hermitian  property  of 
'S'cfdC^mr?  &d,  &rw)  (and  therefore  S^d^mr?  0))  under  experimental  conditions, 
the  estimation  of  Pcfd  is  achieved  by  applying  systematic  phase  corrections 
during  the  computation  of  the  transform  in  Eq.  17.  By  refraining  from  using 
the  signal’s  magnitude  both  the  energy  and  information  content  are  simulta¬ 
neously  preserved  also  permitting  asymmetry  in  Pcfd-  The  starting  point  the 
regular  spin-echo  image:  J^plex((xro,  xpe),  0, /crw)  =  ^"2 ^cfd (^mr ’  ^w)}- 

The  theoretically  real  valued  image  J/Cmr  becomes  complex  valued,  /^omPlex  = 

Ikmr  exp(jZ/^omplex)  due  to  the  systematic  echo  time  shifting  from  the  middle 
of  the  read-out  period  that  creates  a  linear  phase  shift  in  the  physical  read-out 
axis  and  other  effects,  such  as  motion,  along  the  phase-encode  axis.  The  phase 
shifts  are  approximated  by  polynomials  in  each  imaging  axis 

rro(x ro)  ~  Z/CPleX((*ro,  o),  0),  rpe(:rpe)  «  Z/£°™plex((0,  Xpe),  0).  (18) 

The  corrections  are  then  applied  systematically  at  each  value  of  kr>  according 
to: 

hmr(x, kD )  =  7^°“plex(a;,  kD)  exp(-j  (rro(a;ro)  +  rpe(:rPe))).  (19) 
Afterwards,  the  Fourier  transform  in  fo^-space, 

Pcm(x,  wd)  =  ^{hmr(x,  Ad)},  (20) 

is  evaluated  sequentially  in  each  -dimension.  At  each  step  before  executing  the 
transform,  the  Hermitian  property  in  -space  is  re-established,  Ikmr{xi  —  &d)  = 
/^mr(^:^D)  as  demonstrated  in  Fig  2. 

4  Experimental  Setup  and  Analysis  Methods 

A  fixed  baboon  brain  immersed  in  4%  paraformaldehyde  was  used  for  the  exper¬ 
iments.  The  primate  was  prematurely  delivered  on  the  125th  day  and  sacrificed 
on  the  59th  day  after  delivery.  All  animal  husbandry,  handling,  and  procedures 
were  performed  at  the  Southwest  Foundation  for  Biomedical  Research,  San  An¬ 
tonio,  Texas.  Animal  handling  and  ethics  were  approved  to  conform  to  American 
Association  for  Accreditation  of  Laboratory  Animal  Care  (AAALAC)  guidelines. 
Further  details  of  the  preparation  are  explained  in  [11]. 

The  experiments  were  carried  out  on  a  4.7  Tesla  MR  scanner  (Varian  NMR 
Systems,  Palo  Alto,  CA,  USA)  with  a  15  cm  inner  diameter  gradient  system, 
45  Gauss/cm  maximum  gradient  strength  and  0.2  ms  rise  time  using  a  cylindrical 
quadrature  birdcage  coil  (Varian  NMR  Systems,  Palo  Alto,  CA,  USA)  with 
63  mm  inner  diameter. 

CFD-MRI  data  were  obtained  using  the  standard  pulsed-gradient  spin- 
echo  multi-slice  sequence.  The  kmr- space  was  sampled  to  result  in  images  of 
128  x  128  pixels  with  a  field  of  view  64  x  64  mm2  and  0.5  mm  slice  thickness. 
The  &D-space  was  sampled  in  a  uniformly  spaced  Cartesian  grid  in  a  cube 
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Fig.  2.  On  the  top  row  the  Nyquist  plots  of  uncorrected  (left)  and  corrected  (right) 
Scfd  obtained  from  the  experimental  data  described  in  Section  4.  The  plots  show  data 
acquired  at  each  diffusion  gradient  value  km  on  the  complex  plane.  On  the  bottom 
row  are  the  magnitude  and  phase  plots  of  the  data.  Uncorrected  data  (bottom  row  left, 
second  column)  exhibit  a  linear  phase  shift  around  0  frequency,  indicative  of  coherent 
motion.  After  the  phase  corrections  obtained  using  the  polynomial  0.266  km  estimated 
from  the  points  km  =  —6,  0,  6,  12  Gauss/cm,  the  magnitude  is  unchanged  but  signal’s 
imaginary  part  is  smaller  for  the  corrected  values  visible  in  the  vertical  axis  spans  of 
Nyquist  plots  and  the  phase  plots. 

[—30  Gauss/cm,  30  Gauss/cm]3  with  llxllxll  voxels,  i.e.,  6  Gauss/cm  sam¬ 
pling  intervals  at  each  dimension.  In  regular  MRI  terms,  1331  images  per  slice 
were  obtained.  The  repetition  time  Tr  =  Is,  echo  time  Te  =  56.5ms,  diffusion 
pulse  time  offset  A  =  30  ms  and  diffusion  pulse  duration  5  =  15  ms  were  used. 

After  the  acquisition,  the  data  were  transferred  to  a  two  quad  core  2.3GHz 
Intel  Xeon®  cpu  and  8GB  memory  Dell  Precision  Workstation  490  running 
Windows  XP®  64-bit  operating  system.  The  DWI  data  were  placed  in  a  5- 
dimensional  array  in  the  computer  memory  and  the  discrete  Fourier  transform 
was  computed  along  with  the  phase  corrections  as  recounted  in  Section  3.  After¬ 
wards,  isosurfaces  were  computed  as  explained  in  Section  5.  In-house  Mat  lab® 
(Mathworks,  Natick,  MA  USA)  programs  were  used  for  all  of  the  computations 
and  to  display  the  graphics  and  maps. 


5  Visualization  of  the  CFD  Distribution 

In  the  case  of  a  two  dimensional  MR  slice  (Zmr  =  2),  the  domain  of  Pcfd  is  the 
five  dimensional  real  space  Pcfd  :  Mlmr  x  JR3  — »  JR,  which  makes  the  visualiza¬ 
tion  challenging.  The  task  is  realized  by  using  two  components:  the  presenta¬ 
tion  of  the  number  of  magnetic  moments  at  the  initial  position  and  the  three 
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dimensional  distribution  of  the  displacement  integral  values.  The  former  is  ac¬ 
complished  by  using  Pcfd(^,0)  as  the  background  image,  e.g.,  in  Fig.  3,  because 
in  a  restricted  environment  without  coherent  motion,  the  function  Pc$&{x,Wd) 
assumes  its  maximum  at  Wd  =  0  since  the  majority  of  the  magnetic  moments 
do  not  travel  too  far  from  their  initial  position.  If  there  is  coherent  motion,  linear 
portion  of  phase  corrections  rectify  it  anyway. 
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Fig.  3.  The  isosurfaces  (Pcfd  =  0.17)  picked  from  the  area  around  the  corpus  collosum 
(CC)  and  external  capsule  (EC)  junction.  The  area  of  interest  is  indicated  on  the  right 
bottom  corner.  Starting  from  left  bottom  going  clockwise,  the  sample  pixels  of  Fig.  4 
are  from  cerebrospinal  fluid  (CSF),  CC,  White  Matter  (WM)  and  CC  and  EC  junction 
respectively.  The  isosurfaces  for  different  c-values  are  shown  in  Fig.  4. 


In  order  to  present  the  displacement  integral  component  with  the  sole  pur¬ 
pose  of  obtaining  a  single  scale  across  the  image,  Pcfd  is  normalized  at  each 
location  x  :  PcfdOc  Wd)  =  Pcfd(^,  Wd)/Pc fd(x,0).  The  second  part  of  visualiza¬ 
tion  is  achieved  by  overlaying  on  Pcfd(^50)  the  isosurfaces  of  Pcm  with  a  level 
value  c  (0  <  c  <  1),  i.e.  {Wd  E  1R 3  :  PcfdOc  Wd)  =  c}.  For  the  sake  of  an  objec¬ 
tive  assessment,  a  common  c-value  is  adopted  throughout  the  image  after  the 
normalization  creating  a  different  isosurface  at  each  location.  For  biological  tis¬ 
sue,  relevant  information  is  provided  by  the  values  of  Pcfd  away  from  c  =  0,  i.e., 
the  outskirts  of  Pcfd  corresponding  to  the  small  number  of  magnetic  moments 
that  travel  further  away,  portraying  the  microstructure.  The  effect  of  level  value 
choice  is  summarized  in  Fig.  4  displaying  isosurfaces  from  four  different  pixels 
within  the  baboon  brain  marked  in  Fig.  3: 
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1)  Too  high  values  do  not  provide  enough  structural  information  (see  the  first 
rows  in  Fig.  4). 

2)  The  appropriately  informative  value  depends  on  the  properties  of  the  motion 
(thus  of  the  microstructure)  at  a  given  location  (compare  the  columns  of  Fig.  4). 

3)  Too  low  values  force  the  isosurfaces  to  become  extremely  noisy  (see  last  rows 
of  Fig.  4). 
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Fig.  4.  Using  multiple  c-values,  the  isosurfaces  from  the  area  around  the  CC  and  EC 
junction  marked  in  Fig.  3  are  shown  with  different  view  angles  in  JR3  at  each  subfigure 
in  the  quadrants.  Less  and  less  information  is  available  with  an  increasing  c-value.  In 
contrast,  if  the  value  is  chosen  to  be  too  low,  all  of  the  isosurfaces  are  immersed  in 
noise.  Regardless  of  the  c-value,  the  CSF  isosurface  is  noisy  because  the  diffusion  pulse 
times,  which  are  appropriate  for  diffusion  in  tissue,  are  too  long  for  the  displacement 
in  the  fluid. 

In  Figs.  3  and  4,  highly  structured  areas,  e.g.,  the  corpus  callosum  (CC), 
demonstrate  ellipsoidal  shapes.  In  contrast,  at  the  junctions,  the  isosurfaces  in¬ 
dicate  crossing  fibers  (see  the  ‘T5  shaped  isosurface  in  Fig.  4)  at  the  CC-external 
capsule  (EC)  junction.  The  organization  in  white  matter  presents  shapes  in  be¬ 
tween.  These  properties  are  highly  informative  for  tractographic  studies. 

Overall,  the  isosurfaces  are  not  constrained  to  given  forms  like  Gaussians, 
spherical  harmonics  or  to  any  expansions.  In  fact,  they  are  typically  not  even 
symmetric.  They  are  general  and  direct. 
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6  Conclusion  and  Future  Studies 


CFD-MRI  brings  clarifications  to  the  nature  of  the  DW-MR  signal  as  well  as 
conceptual  simplifications  over  existing  methods  [12].  The  narrow  pulse  approxi¬ 
mation  (negligible  duration  of  diffusion  sensitizing  magnetic  field  gradient  pulse) 
is  not  necessary  because  CFD-MRI  uses  the  displacement  integrals,  which  in¬ 
corporate  directly  pulse  durations,  rather  than  the  displacements  themselves.  In 
other  words,  Pc fd  does  not  involve  the  propagator  which  is  the  basis  of  spectral 
methods  [13].  Furthermore,  this  also  removes  the  necessity  to  model  the  motion 
as  a  stochastic  process  relieving  in  conjunction  the  constraints  and  derivations, 
such  as  being  Markovian,  in  the  existing  methods. 

The  high  dimensional  Fourier  relationship  with  the  support  of  Hermitian 
property  preservation  (to  avoid  the  usage  of  signal  magnitude)  results  in  com¬ 
pletely  free  distribution  functions  that  display  asymmetry  as  demonstrated  in 
Figs.  3  and  4  in  an  evidence  based  manner.  Accordingly,  re-establishment  of 
Hermitian  symmetry  requires  future  development  of  more  sophisticated  system¬ 
atic  phase  corrections  starting  from  the  ones  introduced  herein.  From  a  practical 
point  of  view,  this  mathematically  challenging  task  also  requires  a  deeper  investi¬ 
gation  and  better  understanding  of  the  origin  of  DW-MRI  artifacts,  specifically 
the  ones  that  disturb  the  Hermitian  property. 

The  clinical  aim  of  future  studies  is  the  establishment  of  connections  between 
biological  properties  and  the  CFD  distribution  function,  Pc^.  The  challenge  is, 
therefore,  the  classification  of  different  microstructures,  healthy  and  patholog¬ 
ical,  despite  lack  of  analytical  descriptions  of  Pc fd  for  even  simple  geometries. 
A  biological  phantom,  with  known  geometry,  constructed  with  two  nerve  fiber 
bundles  and  agar  gel  has  been  recently  investigated  as  an  initial  step  towards 
that  goal  [14]. 

The  analysis  and  visualization  of  Pc fd  need  to  be  designed  to  define  and 
enhance  CFD-MRI  surrogate  markers  for  diagnosis  and  prognosis  of  disease, 
in  the  mentality  of  the  accomplishments  with  DTI  markers  in  [3].  Naturally, 
these  markers  will  provide  much  deeper  information  than  the  existing  DW- 
MRI  methods.  However,  the  effect  of  integration  on  the  displacement  has  to 
be  carefully  studied,  for  example,  in  order  to  make  axon  caliber  assessments 
for  white  matter  integrity  estimation.  Likewise,  the  observation  that  in  DW- 
MRI  the  displacement  is  from  the  initial  magnetic  moment  position  indicates 
that  fiber  tractography  methods  cannot  exactly  determine  the  location  of  tracts 
within  a  pixel.  It  will  be  necessary  to  introduce  new  tools,  e.g.,  from  differential 
geometry,  for  better  characterization,  categorization  and  visualization  of  three 
dimensional  motion  portion  of  Pcfd- 
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Appendix  3 
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ABSTRACT:  The  addition  of  a  pair  of  magnetic  field  gradient  pulses 
had  initially  enabled  the  measurement  of  spin  motion  to  nuclear  mag¬ 
netic  resonance  (NMR)  experiments.  In  the  adaptation  of  diffusion 
weighted  (DW)-NMR  techniques  to  magnetic  resonance  imaging 
(MRI),  the  taxonomy  of  mathematical  models  is  divided  in  two  catego¬ 
ries:  model  matching  and  spectral  methods.  In  this  review,  the  meth¬ 
ods  are  summarized  starting  from  early  DW  NMR  models  followed  up 
with  their  adaptation  to  DW  MRI.  Finally,  a  newly  introduced  Fourier 
analysis  based  unifying  theory,  so-called  Complete  Fourier  Direct  MRI, 
is  included  to  explain  the  mechanisms  of  existing  methods.  ©  2012 
Wiley  Periodicals,  Inc.  Int  J  Imaging  Syst  Technol,  22,  44-52,  2012;  Published 
online  in  Wiley  Online  Library  (wileyonlinelibrary.com).  DOI  10.1002/ima.22001 

Key  words:  nuclear  magnetic  resonance;  magnetic  resonance 
imaging;  diffusion;  diffusion  weighted  imaging;  Fourier  transform 

I.  INTRODUCTION 

Nuclear  magnetic  resonance  (NMR)  and  investigation  of  diffusion 
phenomenon  have  led  to  important  physical  and  mathematical  dis¬ 
coveries  (Mazo,  et  al.  2002)  in  the  20th  century;  the  discovery  of 
material  properties  via  NMR  spectroscopy  and  advances  in  the 
theory  of  stochastic  processes  to  name  two.  These  research  areas 
come  together  in  the  analysis  of  the  effect  of  spin  diffusion  on 
NMR  signal.  As  the  spins  move  in  a  magnetic  field  gradient  (i.e.,  a 
changing  magnetic  field  strength  dependent  on  position)  their  mag¬ 
netization  vector  changes  individually  according  to  their  paths.  This 
was  accounted  for  in  different  manners  by  the  pioneers  of  NMR: 
Hahn’s  (Hahn  et  al.  1950)  description  of  diffusing  spins’  effects  in 
spin  echo  experiment  followed  by  Carr  and  Purcell’s  (1954)  random 
walk  approach  and  Torrey’s  incorporation  of  magnetization  diffu¬ 
sion  (Torrey,  1956)  to  Bloch  equations.  Torrey  utilized  a  partial  dif¬ 
ferential  equation  (PDE)  framework  instead  of  particle  methods 
used  by  Hahn,  Carr  and  Purcell. 

In  Torrey’s  work  (1956),  the  effects  of  diffusion  on  NMR  signal 
“under  conditions  of  inhomogeneity  in  magnetic  field”  (Torrey, 
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1956)  were  mathematically  formulated.  Shortly  thereafter,  moving 
one  step  further  Stejskal  and  Tanner  intentionally  introduced  a  pair 
of  (homogenous)  magnetic  field  gradients  into  the  spin  echo  experi¬ 
ment  with  the  purpose  of  accurately  measuring  the  scalar  diffusion 
coefficient  of  the  sample  under  investigation  (Stejskal  and  Tanner, 
1965).  When  spins  diffuse  in  a  magnetic  field  gradient  their  trans¬ 
verse  magnetizations  disperse  because  the  gradients  or  inhomoge¬ 
neity  creates  a  different  magnetic  field  strength  at  each  location. 
This  changes  the  rotational  speed  of  each  magnetization  vector 
according  to  its  spin’s  diffusion  path.  Each  spin  possess  a  different 
path  thereby  creating  the  dispersion  of  the  magnetization.  In  the 
particle  method  framework,  this  principle  was  also  analyzed  by 
McCall  et  al.  (1963)  with  a  statistical  approach  that  refined  Hahn’s 
(1950)  magnetization’s  phase  distribution  formulation  but  departing 
from  the  random  walk  model  of  Carr  and  Purcell  (1954). 

Although  the  pair  of  magnetic  field  gradient  pulses  are  conven¬ 
tionally  named  as  diffusion  sensitizing  gradients  in  the  pulsed  gradi¬ 
ent  spin  echo  experiment  (PGSE)  of  Stejskal  and  Tanner,  they  do  in 
fact  encode  spin  motion  regardless  of  its  coherence  or  lack  thereof 
(Nalcioglu  and  Cho,  1987).  In  a  sense,  diffusion  NMR  has  introduced 
the  motion  information  as  an  additional  dimension  before  magnetic 
resonance  imaging  (MRI)  increased  dimensionality  with  addition  of 
position  of  spins  as  an  independent  variable.  Nevertheless,  the  tech¬ 
niques  involving  diffusion  sensitizing  gradients  whether  for  NMR 
spectroscopy  or  MRI  are  called  diffusion  weighted  (DW)  methods. 

Decades  later,  the  DW  technique  has  been  adapted  to  MRI  with 
great  success,  with  early  detection  of  ischemia  (Moseley  et  al.,  1990) 
being  a  very  important  example.  Investigations  were  based  on  the  esti¬ 
mation  of  a  scalar  descriptor,  the  so-called  apparent  diffusion  coeffi¬ 
cient  (ADC),  which  represents  the  equivalent  of  isotropic  diffusivity 
in  an  anisotropic  environment.  In  other  words,  environments  allowing 
unrestricted  and  restricted  motion  of  molecules  are  compared  with 
ADC.  DW  magnetic  resonance  (MR)  measurement  of  anisotropy 
caused  by  three-dimensional  microstructure’s  alteration  of  spin 
motion  raised  an  opportunity  to  infer  more  detailed  information,  spe¬ 
cifically  about  biological  tissue.  Furthermore,  researchers  took  advant¬ 
age  of  MRI  hardware’s  capability  of  producing  vectorial  magnetic 
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field  gradients.  They  developed  models  to  express  characteristics  of 
spin  motion  as  three  dimensional  constructs  which  traversed  beyond 
the  scalar  (apparent)  diffusion  coefficient’s  potential.  The  most  widely 
known  of  these  models  are  perhaps  the  q-space  formalism  (Callaghan, 
1991)  and  diffusion  tensor  imaging  (DTI)  (Basser  et  al.,  1994;  Mat- 
iello  et  al.,  1994).  Whereas  the  former  builds  a  Fourier  framework 
between  DW-MR  signal  and  spin  displacement  properties,  the  latter 
projects  the  signal  onto  an  ellipsoid,  mathematically  described  by  a 
symmetric  tensor  of  rank  two  (which  is  a  quadratic  form). 

These  developments  have  pushed  the  initial  utilization  of  ADC 
to  different  crucial  areas  in  research  and  clinical  imaging:  early 
clinical  diagnosis  of  ischemia  (Baird  and  Warach,  1998),  cancer  di¬ 
agnosis  (Song  et  al.,  2002;  Xu  et  al,  2009)  follow-up  on  treatment, 
preoperative  and  postoperative  assessment  for  different  organs 
[e.g.,  fiber  tracking  (Conturo  et  al.,  1999;  Mori  and  van  Zijl,  2002) 
before  brain  surgery  (Warfield,  2003)],  monitoring  of  neurological 
diseases  and  disorders  (Budde  et  al.,  2008;  Ciccarelli  et  al.,  2008), 
neonatal  development  (McKinstry  et  al.,  2002)  and  traumatic  brain 
injury  (Mac  Donald  et  al.,  2011).  In  consequence,  DW-MRI  is  now 
an  indispensable  and  versatile  tool,  widely  used  in  research  with  an 
increasing  number  of  accepted  clinical  applications. 


II.  DIFFUSION  WEIGHTED  MODELS  FOR 
NMR  SPECTROSCOPY 

A.  PDEs.  The  mathematical  treatment  of  the  diffusion  MR  signal 
is  confined  traditionally  to  a  PDE  framework,  which  is  used  for 
both  diffusion  (Crank,  1975)  and  Bloch  equations.  The  starting 
point  is  Fick’s  first  equation  that  describes  the  rate  of  transfer,  /,  of 
diffusing  magnetization,  M,  through  the  unit  area  in  each  direction 
where  D  is  the  diffusion  coefficient: 


/  =  -DVM. 


(i) 


Whereas  Torrey  (1956)  leaves  the  scalar  diffusion  coefficient,  D,  as 
position  dependent  (thereby  covering  a  more  general  case),  Stejskal 
and  Tanner  (1965)  following  Abragam’s  formulation  (Abragam,  1961) 
treat  it  as  position  independent  in  Fick’s  second  equation  (Crank, 
1995)  allowing  the  commutation  of  operators  in  the  diffusion  equation, 


^  =  V.DVM  ^^  =  DV2M 
at  at 


(2) 


( t  denotes  time).  This  formulation,  valid  for  unrestricted  diffusion, 
as  in  a  liquid  sample,  describes  the  Brownian  motion.  Accordingly, 
Brownian  motion  is  mathematically  modeled  by  Wiener  process, 
which  is  defined  as  the  stochastic  process  with  Gaussian  distributed 
independent  increments.  By  this  property,  the  Wiener  process 
equivalently  satisfies  the  Markovian  property  (Shiryaev,  1996).  In 
other  words,  the  spins  in  Brownian  motion  are  memoryless,  their 
future  positions  do  not  depend  on  their  position  history.  On  the 
other  hand,  the  obstacles  posed  by  microstructure  determine  the 
future  of  spin  motion  according  to  past  spin  interactions  with  micro¬ 
structure.  Therefore,  the  inclusion  of  Markovian  property  is  not  nec¬ 
essarily  adequate  in  describing  molecular  motion  in  a  restricted 
environment  where  past  positions  shape  future  motion. 


In  the  investigations  of  microstructured  environments,  this  is  the  reason  the  term 
‘apparent’  appears  in  front  of  diffusion  coefficient,  which  is  defined  for  environments 
that  allow  unrestricted  or  isotropic  motion. 


s. 


Single  Point  Measurement 


Figure  1.  Spin  echo  (top)  and  pulsed  gradient  spin  echo  experiments. 
Diffusion  sensitizing  gradients  are  placed  around  the  180  degree  radio  fre¬ 
quency  (RF)  pulse.  The  duration  and  the  separation  of  the  pulses  are  8  and 
A  respectively.  The  echo  time  is  denoted  by  Te  and  the  measurement  is 
the  attenuation  of  the  single  point,  the  maximum  of  FID.  [Color  figure  can 
be  viewed  in  the  online  issue,  which  is  available  at  wileyonlinelibrary.com.] 


Nevertheless,  the  Bloch-Torrey  equation,  which  phenomenologically 
describes  NMR,  takes  this  form  after  neglecting  relaxation  effects: 

^  =  yMXB  +  DV2M ,  (3) 

with  y  and  B  denoting  the  gyromagnetic  ratio  and  the  magnetic  field 
respectively.  Proceeding  with  the  solution  of  Bloch-Torrey  equation 
for  PGSE,  Stejskal  and  Tanner  rigorously  obtained  the  characteriza¬ 
tion  of  the  DW  signal  attenuation  of  the  free  induction  decay’s 
(FID)  magnitude  as: 


Si=S0Gv(-bD).  (4) 

Here  St  denotes  the  magnitude  of  the  FID  corresponding  to  the  acqui¬ 
sition  using  the  ith  diffusion  sensitizing  gradient  magnitude  with  i  = 
0  indicating  the  absence  of  diffusion  sensitizing  gradients  (see  Fig.  1) 
and  b  is  derived  using  the  duration,  A,  the  separation,  A,  times  of  the 
(rectangular)  diffusion  sensitizing  gradients  and  the  magnitude  of  the 
diffusion  sensitizing  gradient  pulse  g  (Stejskal  and  Tanner,  1965): 

6  =  7282  (A-lsy.  (5) 

Basically,  the  attenuation  is  characterized  as  a  Gaussian  function  of 
diffusion  gradient  strength  g.  The  diffusion  coefficient  and  the  dif¬ 
fusion  gradient  times  define  the  shape  of  the  Gaussian. 


B.  Displacement  Probabilities.  Before  the  appearance  of  Tor- 
rey’s  PDE  approach  (Torrey,  1956),  initial  modeling  of  diffusion 
NMR  in  spin  echo  experiments  by  Hahn  (1950)  was  based  on  the 
expression  of  the  distribution  of  individual  spins’  transverse  magnet¬ 
ization  phase.  Carr  and  Purcell  (1954)  interpreted  diffusion  using  a 
random  walk  model  to  arrive  at  Hahn’s  result  for  signal  attenuation: 

S;  =  Soexp(-^Y^o)-  (6) 

Unlike  in  the  case  of  PGSE  experiment,  g  represents  here  a  constant 
magnetic  field  gradient  rather  than  the  strength  of  pulsed  gradients 
of  Figure  1. 
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Around  the  publication  time  of  Stejskal  and  Tanner’s  seminal 
work  (Stejskal  and  Tanner,  1965),  another  less  noticed  work  treat¬ 
ing  diffusion  NMR  signal  via  particle  methods  was  published  by 
McCall  et  al.  (1963).  The  work  departed  from  the  random  walk  for¬ 
mulation  of  Carr  and  Purcell  (1954),  which  records  the  phase  of  the 
transverse  magnetization  during  each  random  walk  step  in  a  mag¬ 
netic  field  gradient.  The  analysis  of  McCall  et  al.  (1963)  returned  to 
the  initial  statistical  treatment  of  Hahn  (1950)  describing  the  mag¬ 
netization  phase  distribution  caused  by  spin  motion.  Despite  minor 
mathematical  discrepancies^  the  fruitful  idea  refines  Hahn’s  (1950) 
analysis  by  expressing  the  individual  phase  directly  as  a  function  of 
spin  displacements  and  obtains  exactly  Eq.  (6). 

More  than  two  decades  after  the  publication  of  McCall  et  al. 
(1963),  Callaghan  modeled  the  PGSE  signal  using  techniques  of 
scattering  theory  in  the  spirit  of  particle  methods.  In  the  model,  an 
analogy  between  the  probability  of  “finding  any  scattering  center” 
(Callaghan,  1991)  and  the  probability  of  spin  displacements  creates 
the  so-called  q-space  model.  With  the  narrow  pulse  assumption, 
i.e.,  negligible  diffusion  gradient  pulse  duration  in  comparison  with 
the  separation  time,  8  <<  A,  Callaghan  computes  the  amplitude 
of  the  echo  at  Te  as 

E(Gd)  =  J  p(x)  J  Ps(d,  A|x)exp(/y8G£>'(x/  —  x))dx'dx  (7) 

with  p  and  GD  denoting  the  spin  density  and  the  three  dimensional 
diffusion  sensitizing  magnetic  field  gradient  vector  respectively.  Ps 
is  the  probability  of  a  particle  initially  at  x  to  be  found  at  x'  after  the 
pulse  separation  time,  A.  Relying  on  Markovian  property,  Calla¬ 
ghan  rewrites  Eq.  (7)  using  the  “dynamic  displacement”:  R  =  x'  — 
x,  and  the  average  propagator  (Callaghan,  1991),  P(R,A),  which 
“gives  the  average  probability  for  any  particle  to  have  a  dynamic 
displacement  R  over  a  time  A”  (Callaghan,  1991).  Accordingly  the 
signal  at  echo  time  is  calculated  as: 

E(q)  =  J  P(R,  A)exp(/2jt  q  •  R)dR,  (8) 

where  q  =  (27r)_1y8G,  giving  the  name  q-space. 

The  derivations  of  these  results  are  strongly  dependent  on  the 
single  point  measurements  at  echo  time  (see  Fig.  1).  However,  the 
same  reasoning  is  not  applicable  to  DW-MRI  for  the  reasons  that 
will  be  explained  in  Section  4. 

Finally,  it  is  worth  pointing  out  that  particle  methods  and  PDE’s 
are  equivalent  mathematically  (Feller,  1967).  Consequently,  both 
PDE  and  particle  methods  based  derivations  are  valid  for  DW 
methods. 

III.  DW  MRI  MODELS 

In  comparison  with  NMR  spectroscopy,  MRI  scanner  hardware  cre¬ 
ates  three  dimensional  magnetic  field  gradients  that  allow  the  mea¬ 
surement  of  spin  density  as  a  function  of  position.  The  spin  density, 
usually  weighted  by  relaxation  times  (7\,  T2)  is  presented  as  an 
image.  In  other  words,  as  an  improvement  to  obtaining  an  average 
measurement  from  the  whole  sample  of  NMR  spectroscopy,  MRI 
augments  the  signal  dimensionality  by  adding  positional  informa- 

'The  assertion  that  the  integrals  of  sample  paths  of  a  Gaussian  distributed  stochas¬ 
tic  process  are  Gaussian  does  not  exist  in  the  reference  provided  in  the  manuscript 
(Munroe,  1951.).  In  fact,  its  proof  does  not  exist  in  the  literature  to  the  best  of  the 
authors’  knowledge. 


tion.  Beginning  in  the  late  1980s,  diffusion  methods  were  adapted 
to  MRI  using  the  mathematical  models  of  DW-NMR  described  in 
Sections  II. A.  and  II.B.  An  important  leap  occurred  when  research¬ 
ers  moved  from  the  investigations  of  the  scalar  quantity  (ADC)  to 
the  estimation  of  more  general  three  dimensional  geometric  con¬ 
structs.  For  example,  DTI  estimates  ellipsoids  defined  by  a  quad¬ 
ratic  form  (diffusion  tensor),  D,  which  is  basically  a  symmetric  3  X 
3  matrix.  The  scalar  quantity  ADC  is  equal  to  the  mean  of  the 
eigenvalues  of  the  matrix  D. 

DW-MRI  methods  can  be  grouped  under  two  general  categories: 

1.  Model  Matching  Methods  initiated  by  DTI  (Basser  et  al., 
1994;  Matiello  et  al.,  1994)  and  refined  with  high  angular  re¬ 
solution  diffusion-weighted  imaging  (HARDI)  (Frank, 
2001),  composite  hindered  and  restricted  model  of  diffusion 
(CHARMED)  (Assaf  and  Basser,  2005),  spherical  deconvo¬ 
lution  (Tournier  et  al.,  2004),  diffusion  orientation  transform 
(DOT)  (Ozarslan  et  al.,  2006)  two  versions  of  the  generalized 
DTI  (GDTI)  (Ozarslan  and  Mareci,  2003;  Liu  et  al.,  2004) 
and  diffusional  kurtosis  imaging  (DKI)  (Jensen  et  al.,  2005). 

2.  Spectral  methods  originating  from  Callaghan’s  q-space 
(Callaghan,  1991)  followed  by  the  diffusion  spectrum  imag¬ 
ing  (DSI)  (Wedeen  et  al.,  2005),  Q-ball  imaging  (QBI)  (Tuch 
et  al.,  2003). 


A.  Model  Matching  Methods.  DW-MRI  model  matching  meth¬ 
ods  are  direct  adaptations  of  signal  attenuation  modeling  in  DW- 
NMR  spectroscopy.  The  models  assume  that  images  obtained  with 
diffusion  sensitizing  gradients  demonstrate  different  attenuation 
levels  at  different  pixels  as  a  function  of  diffusion  gradient  vectors. 
This  reflects  characteristics  of  spin  motion  dictated  by  the  micro¬ 
structure.  The  model  matching  method  is  described  by  the  follow¬ 
ing  generic  equation  adopted  from  Eq.  (4): 

Ii  =  /o  exp (— //(GdO).  (9) 

Here,  the  model  is  no  longer  for  the  FID  but  its  Fourier  transform, 
which  is  in  the  image  domain:  7/  is  the  intensity  (magnitude  in  most 
cases)  of  a  given  pixel  from  the  image  obtained  with  the  ith  diffu¬ 
sion  gradient  vector  GD/.  The  methods  are  defined  by  the  functional 
presentation  of  the  model,  77.  The  data  are  projected  onto  77,  which 
also  incorporates  the  unknown  to  be  estimated,  namely  the  spin 
motion  descriptor. 

A.l.  DTI.  In  the  original  formulation  of  model  matching  pio¬ 
neering  work,  DTI  (Basser  et  al.,  1994;  Matiello  et  al.,  1994),  the 
function  is 

77  =  7?(Gd/,  6,  A)  :D.  (10) 

The  descriptor  is  the  diffusion  tensor  D  to  be  estimated  from  the  (now 
3  dimensional)  Bloch-Torrey  equation  [Eq.  (3)].  B  is  a  matrix  valued 
function  of  GD/,  8,  A  (Matiello  et  al.,  1994)  and  it  multiplies  D  ele- 
ment-by-element.  This  formulation  is  refined  to  a  natural  linear  alge¬ 
braic  framework  by  Papadakis  et  al.  (1999)  who  noticed  that  in  fact 

B(GDh  8,  A)  :  D  =  y 2b,GTDiDGm.  (11) 

The  constant  coefficient  b,  =  82  (A  —  1 8) ,  introduced  later  in 
(Ozcan,  2005)  delineates  the  gradient  duration  and  separation  times 
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Figure  2.  DW  images  from  a  PGSE  experiment  of  an  ex-vivo  baboon  brain  (Ozcan  201 0;  Ozcan,  201 1)  shown  on  top  of  the  diffusion  sensitizing 
gradient  vectors  (orange  vectors  with  blue  ends)  of  the  DGS.  The  image  at  the  origin  is  acquired  without  diffusion  gradients  and  therefore  is  unat¬ 
tenuated.  The  attenuation  is  highest  on  the  fibers  parallel  to  the  diffusion  gradient  vector  as  shown  in  two  zoomed  images  on  the  left  and  right 
bottom.  [Color  figure  can  be  viewed  in  the  online  issue,  which  is  available  at  wileyonlinelibrary.com.] 


(8,  A)  from  their  vector  part  GD/  in  contrast  to  b  of  Eq.  (5).  Conven¬ 
tional  diffusion  gradient  schemes  (DGS,  the  list  of  diffusion  sensi¬ 
tizing  gradient  vectors  used  in  a  DW  experiment)  are  specified 
using  a  single  lvalue  for  the  entire  diffusion  gradient  vectors  GD/ 
thereby  confining  all  of  them  onto  a  sphere.  The  introduction  of  bt 
made  possible  the  design  of  new  DGS  via  optimization  methods.  In 
(Ozcan,  2011),  the  effect  of  imaging  gradients  is  minimized  with 
DGSs  containing  vectors  of  different  magnitudes  optimally  placed 
with  respect  to  the  imaging  magnetic  held  gradients. 

At  a  basic  level,  Eq.  (11)  implies  that  D  estimated  from  DTI 
experiments  is  forcefully  a  symmetric  matrix  because  for  any  vector 
G,  GtDG  =  GtDtG  and 

2 GtDG  =  Gr(D  +  Dt)G  +  Gt{D  -  Dt)G  =  Gt(D  +  Dt)G. 

(12) 


In  other  words,  DTI  can  only  measure  the  symmetric  part  of  the  dif¬ 
fusion  tensor:  \  (D  +  Dr).  Symmetric  matrices  form  a  subspace  and 
3X3  symmetric  matrices  are  represented  as  six  dimensional  vec¬ 
tors,  d  =  [di,  d2 ,  d3,  d4,  d5,  d6\ : 


D  = 


d\ 

d4 

d6 


d4  dfo 

d2  d$ 
d$  d3 


(13) 


This  allows  rewriting  Eq.  (11)  in  a  linear  form  following  the  steps 
of  Papadakis  et  al.  (1999)  and  the  refinements  in  (Ozcan,  2005): 


72/j,Gq,DGd,  =  y2b,Vid  (14) 


with 


Vi  —  [§ixi  §iyi  Sizi^SixSiyi^SiySizi^gixg^  0^) 


where  Gd,  =  [gix,giy,giz]- 

To  solve  for  six  elements  of  d,  a  minimum  of  six  DW  images  plus 
one  unattenuated  image,  7o,  are  acquired,  m  >  6  (see  Fig.  2)  and  the 
logarithm  of  Eq.  (9)  is  taken  to  give  a  set  of  linear  equations: 


vi 

~  ln(/0)  ln(/i) 

y  2btVgd  =  y 2bt 

d  = 

_vm_ 

_ln(/0)  —  ln(/m)  _ 

The  list  of  the  gradient  vectors  used  for  the  acquisition  of  DW 
images  is  called  the  DGS  (Ozcan,  2010a,  2010b,  2011).  An  example 
of  DGS  with  the  corresponding  images  is  shown  in  Figure  2.  It  fol¬ 
lows  from  elementary  linear  algebra  that  for  Eq.  (16)  to  have  a 


unique  solution,  the  coefficient  matrix,  Vg 


vi 


must  have  full 


_vm_ 

rank.  Vg  is  a  function  of  DGS  vectors  [see  Eq.  (15)]  and  therefore 
the  vectors  must  be  chosen  appropriately  to  guarantee  the  unique¬ 
ness  of  the  solution.  A  set  of  mathematical  necessary  conditions  for 
the  choice  of  diffusion  gradients  that  would  guarantee  the  full  rank 
condition  were  given  in  (Ozcan,  2005)  a  decade  after  the  introduc¬ 
tion  of  DTI. 


Once  the  solution  for  d  is  found,  the  symmetric  matrix  D  is 
reconstructed  as  in  Eq.  (13)  at  each  pixel.  A  standard  theorem  in 
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linear  algebra  proves  that  symmetric  matrices,  thus  the  diffusion 
tensor,  have  orthogonal  eigenvectors  and  real  eigenvalues  (ki  > 
7^2  >  A,3).  The  setup  is  general  in  the  sense  that  ADC  can  be  calcu¬ 
lated  as  the  mean  of  the  eigenvalues. 

In  this  three-dimensional  structure,  whereas  the  eigenvectors 
indicate  the  principal  directions  of  the  motion,  corresponding  eigen¬ 
values  describe  the  tendency  of  the  spins  to  move  in  each  direction. 
If  all  the  eigenvalues  are  equal,  the  medium  is  isotropic  as  in  a  liq¬ 
uid.  However,  if  the  micro  structure  favors  motion  in  a  given  direc¬ 
tion,  causing  anisotropy,  the  eigenvalues  will  reflect  that  by  differ¬ 
ing  from  each  other.  Functions  of  the  three  eigenvalues  describing 
their  dispersion  are  called  anisotropy  measures,  fractional  anisot¬ 
ropy  being  one  of  the  most  used  ones  (Basser  and  Pierpaoli,  1996): 

\!  —  X2)2  +  (^2  —  ^3)2  +  (^3  —  ^l)2 

FA  =  A -  - -.  (17) 

Y  2(^1  +  ^2  +  ^3) 


Usually,  an  anisotropy  map  describing  the  microstructure  is  shown  as 
an  image  with  the  value  of  anisotropy  index  at  each  pixel.  The  direc¬ 
tional  information  might  be  added  by  displaying  the  eigenvector  cor¬ 
responding  to  the  largest  eigenvalue  at  each  pixel  as  a  “whisker" 
(thus  the  name  whisker  plot).  The  directionality  is  also  shown  with  a 
color  image.  A  different  color  channel  is  assigned  to  each  direction  to 
display  three  components  of  the  eigenvector  accordingly.  There  are 
different  applicable  schemes  depending  on  the  properties  of  the  spin 
motion  (Pajevic  and  Pierpaoli,  1999).  The  combined  directional  and 
motion  tendency  can  also  be  shown  by  overlaying  ellipsoids  defined 
by  the  diffusion  tensor  (xTDx  =  constant)  on  each  pixel  of  an  ana¬ 
tomic  image.  The  equation  of  an  ellipse  requires  that  the  eigenvalues 
of  D  are  all  greater  than  or  equal  to  0.  Otherwise,  the  diffusion  equa¬ 
tion,  Eq.  (2),  contradicts  the  conservation  of  mass.  This  physical  con¬ 
straint  is  sometimes  violated  in  the  estimation  obtained  from  experi¬ 
mental  DTI  data.  The  reasons  for  the  occurrences  of  negative  eigen¬ 
values  from  experimental  data  constitute  currently  an  open  question. 


B.  High  Order  Model  Matching  Methods.  DTI  can  be  consid¬ 
ered  as  a  second  order  approximation  to  the  attenuation  function  H  of 
Eq.  (9).  Different  expansions  of  the  attenuation  function  result  in  dif¬ 
ferent  model  matching  methods.  In  HARDI  (Frank,  2001),  the  attenu¬ 
ation  is  expanded  using  the  spherical  harmonics  Yf  (Frank,  2001) 

00  1 

w(GDl(e,  <p))  =  Y2  Y2  a"»Fr(e.  <p)  (is) 

1=0  m——  1 


In  DKI  (Jensen  et  al.,  2005),  a  higher  order  expansion  using  the 
concept  of  excess  kurtosis  is  implemented.  H  is  expanded  by  adding 
the  kurtosis  tensor  to  the  diffusion  tensor  D  (Jensen  et  al.,  2005): 

ADC 

H(Gm)  =  ^btG^DGui - 7—  (Go/ ,  Gx>/ ,  Gx>/ ,  Go/)  • 

o 

(19) 

The  kurtosis  tensor  ZJ(4)  is  a  symmetric  tensor  of  rank  4  and  the 
method  is  a  generalization  of  DTL*  The  expansion  strategy  of  Eq. 
(19)  appeared  before  the  work  of  Jensen  et  al.  (2005).  Ozarslan  and 
Mareci  (2003)  integrated  higher  rank  tensors  into  the  Bloch-Torrey 
equations.  The  odd  ranked  tensor  elements  result  in  negative  diffu¬ 
sion  coefficients  (Ozarslan  and  Mareci,  2003).  Ozarslan  and  Mareci 
(2003)  argue  that  “since  negative  diffusion  coefficients  are  non¬ 
physical,  the  rank  is  forced  to  be  an  even  number”  (see  also  the  dis¬ 
cussion  about  the  negative  eigenvalues  in  Section  III.A.l.).  Conse¬ 
quently,  the  generalized  diffusion  tensor  (GDT)  approximation  of  H 
is  given  as: 


H(GDi)=y2btJ2D(2n) >  (2°) 

n=  1 


The  n  arguments  each  of  D(n)  s  operates  on  is  GD/,  e.g.,  in  the  last  term 
of  Eq.  (19)  where  n  =  4.  About  a  year  after  the  publication  of  (Ozar- 
slan  and  Mareci,  2003),  Liu  et  al.  (2004)  published  another  GDT  imag¬ 
ing  model  that  incorporates  higher  order  tensors  of  all  ranks: 


00 

H(GDi)  = 

n= 2 


(21) 


where  j  =  and  b\n^  ’s  are  higher  order  integrals  of  the  diffusion 
sensitizing  gradient’s  time  course.  The  inclusion  of  higher  rank  ten¬ 
sors  is  made  possible  by  allowing  the  signal  to  be  complex  valued. 
In  Eq.  (19),  odd  ranked  tensors  form  the  imaginary  part  of  the  sig¬ 
nal.  In  addition,  they  are  asymmetric  terms  and  this  covers  physi¬ 
cally  a  more  general  situation. 

The  general  aim  of  these  model  matching  methods  is  to  estimate 
the  descriptors  of  diffusion  that  would  minimize  the  model  match¬ 
ing  error  originating  from  the  measurement  values,  /1,  ...,/m: 


-/o 


exp  (— //(Gm))' 

exp(— //(GDm)) 


(22) 


where  (0,  (p)  are  the  azimuthal  and  polar  angles  in  the  spherical 
coordinates.  The  coefficients,  aim,  are  computed  from  Eq.  (18)  using 
the  experimental  data.  The  odd  numbered  harmonics  /  =  1,  3,  5. . . 
are  discarded  on  the  basis  that  they  are  asymmetric  (Frank,  2001). 

A  generalization  of  HARDI  is  provided  by  Tournier  et  al. 
(2004)  in  spherical  deconvolution,  utilizing  concepts  from  linear 
systems  theory  with  the  use  of  a  response  function  (also  known  as 
point  spread  function  or  transfer  function  in  control  theory)  in 
spherical  coordinates.  The  response  function,  equivalent  to  impulse 
response  of  a  linear  time  invariant  dynamical  system,  describes  the 
signal  originating  from  a  single  nerve  fiber.  To  assess  the  distribu¬ 
tion  of  multiple  fibers  within  a  voxel,  the  full  set  of  measurements 
is  deconvolved  with  the  response  function.  The  result  is  the  fiber 
orientation  density  function. 


The  solution  of  the  minimization  problem  of  Eq.  (22)  for  each 
model  matching  method  results  in  the  determination  of  a  number  of 
unknowns  (e.g.,  6  for  DTI)  specified  by  the  model  equations  pre¬ 
sented  in  this  section.  Naturally  this  number  is  also  the  minimum 
number  of  diffusion  gradient  vectors  in  DGS.  Furthermore,  in  a  par¬ 
allel  observation  to  the  DTI  estimation  equations,  the  conditions  on 
DGS  guaranteeing  the  uniqueness  of  solution  for  different  model 
matching  methods  have  not  been  investigated  and  remain  as  an 
open  and  untouched  subject  of  discussion. 

Tor  example,  D  is  a  rank  two  tensor  and  both  of  its  arguments  are  equal  to  GD;. 
Therefore,  it  is  a  quadratic  form  that  can  be  represented  as  a  matrix.  Its  value  at  GD,  is 
found  by  left  and  right  vector-matrix  multiplication  i.e.,  D(Gd„  GD;)  =  GDiTDGD;. 
Higher  rank  tensors  have  more  complicated  representations  than  vector-matrix 
multiplication. 
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In  addition  to  the  aforementioned  methods  that  are  using  the 
attenuation  model  of  Eq.  (9),  there  are  two  techniques  that  are  com¬ 
partmentalizing  and  approximating  the  q-space  signal  of  Section 
II.B.  The  first  one,  composite  hindered  and  restricted  model  of  dif¬ 
fusion  (CHARMED)  (Assaf  and  Basser,  2005),  divides  the  signal 
into  two  compartments  as  “hindered  water  diffusion  in  the  extra- 
axonal  space  and  restricted  water  diffusion  in  the  intra-axonal 
space.”  Accordingly  the  q-space  signal  originating  from  multiple 
compartments  at  the  echo  time  is  given  as 

E{q)  =  YfA  +  YS&  (23) 

i=  1  7=1 

where  4  4L  is  the  number  of  distinct  hindered  compartments  and  N  is 
the  number  of  distinct  restricted  components”  (Assaf  and  Basser, 
2005),  /h  and  fj.  are  volume  fractions  of  the  compartments.  Each  of 
the  compartments  are  further  dissected  to  divide  the  signal  into  par¬ 
allel  and  perpendicular  (e.g.,  to  axons)  sections  to  obtain  a  final 
model.  The  second  technique,  DOT  of  Ozarslan  et  al.  (2006),  takes 
the  different  path  of  expanding  theoretically  the  Fourier  Transform 
in  Eq.  (9)  using  spherical  harmonics.  The  propagator,  P  of  Eq.  (8), 
is  expressed  in  the  basis  formed  by  the  spherical  harmonics.  The 
expression,  in  essence,  is  the  same  as  Eq.  (18)  but  it  is  the  propaga¬ 
tor,  rather  than  the  attenuation,  that  is  approximated. 

C.  Spectral  Methods.  The  basis  of  the  spectral  methods  for 
DW-MRI  is  the  adaptation  of  NMR  q-space  equation,  Eq.  (8),  to 
MRI.  (Callaghan,  1991): 

S(k,q)  =  J  p(x)E(q)  exp(z'27t  k  ■  x)dx,  (24) 

where  k  is  the  usual  MRI  k- space  variable  and  E{q)  is  the  adapta¬ 
tion  of  its  NMR  definition  in  Eq.  (7) 

E(q)  =  Jps{x'A\X)zMn*GD  ■(x/—x))dxr.  (25) 

Clearly,  in  this  equation,  E(q)  is  a  function  of  position  x.  Accord¬ 
ingly,  it  must  be  taken  into  account  for  the  evaluation  of  the  integral 
in  Eq.  (24).  To  obtain  Eq.  (24),  on  p.  440  of  (Callaghan,  1991)  it  is 
asserted  that  “despite  the  fact  that  both  Ps  and  p  may  depend  on  x, 
the  effect  of  the  PGSE  sequence  is  quite  separable. . .”.  It  is  also 
indicated  that  the  signal  comes  from  the  finite  volume  element  cor¬ 
responding  to  the  imaging  voxel.  In  consequence,  E(q)  should  be 
seen  as  “the  Fourier  transform  of  local  dynamic  profile”  according 
to  Eq.  (8.7),  p.  440,  (Callaghan,  1991): 

E(q)  =  J  P(R ,  A)  exp(/27t  q  ■  R)dR.  (26) 

The  goal  of  the  method  is  the  reconstruction  of  the  average  propa¬ 
gator  using  the  Fourier  transform  from  the  measurements  after  their 
normalization.  This  is  explicitly  formulized  in  diffusion  spectrum 
imaging  (DSI)  method  of  Wedeen  [Eq.  (4),  (Wedeen  et  al.,  2005]: 

P{R,  A)  =  /„  : l{2ny3  J  |/(<7)|exp(— 2/jr  q  .  R)dq.  (27) 

Here,  I(q)  is  the  (complex)  pixel  intensity  corresponding  to  diffu¬ 
sion  sensitizing  gradient  that  produces  q.  In  other  words,  I(q)  =  7/ 


which  is  the  image  obtained  with  GD/  such  that  q  =  (27r)_1y8GD/. 
The  magnitude  is  taken  “to  exclude  phase  shifts  arising  from  tissue 
motion”  (paragraph  preceding  Eq.  (4),  (Wedeen  et  al.,  2005]).In 
DSI  the  main  interest  is  the  angular  structure  of  the  diffusion  spec¬ 
trum.  The  final  product  is  the  orientation  distribution  function 
defined  as  ’a  weighted  radial  summation  of  P(R,  A)  (Wedeen  et  al., 
2005): 

ODF(u)  =  J  P(ru,A)r2dr  with  ||«||  =  1  for  r  >  0.  (28) 

It  is  important  to  note  the  derivations  in  (Wedeen  et  al.,  2005)  rely 
on  Markovian  property  and  the  narrow  pulse  approximation.  More¬ 
over,  the  magnitude  of  the  signal  is  used  on  the  basis  that  “the  MRI 
signal  is  positive  for  any  type  of  spin  motion  without  net  flux” 
(Wedeen  et  al.,  2005).  Under  these  assumptions,  the  model  is  not 
completely  unconstrained. 

The  Fourier  reconstruction  in  general  requires  a  large  number  of 
DW  acquisitions  compared  with  model  matching  methods.  In  DSI, 
the  sampling  of  q-space  is  realized  within  a  sphere.  The  issue  of 
large  number  of  acquisition  is  addressed  by  the  use  of  the  so-called 
QBI  technique  (Tuch  et  al.,  2003).  As  in  DSI,  in  QBI  the  quantity 
of  interest  is  the  orientation  distribution  function  which  is  defined 
slightly  different  that  DSI-ODF: 

v| j(u)  =  J  P(ru,A)dr  with  ||w||  =  1  for  r  >  0.  (29) 

By  utilizing  of  the  Funk-Radon  transform,  QBI  estimates  \| /  with 
measurements  obtained  on  a  ball  (hence  the  method’s  name,  QBI) 
in  q-space,  with  significantly  smaller  number  of  acquisitions. 

IV.  A  HIGHER  DIMENSIONAL  FOURIER  RELATIONSHIP 

The  methods  described  in  previous  sections  are  just  representatives 
of  many  existing  variations  of  DW-MR  methods  and  techniques  that 
were  not  listed  here  due  to  space  restrictions.  An  exhaustive  list  can 
be  found  in  the  publication  by  Jian  et  al.  (2007).  The  sheer  number 
of  methods  could  be  considered  as  an  indication  of  dissatisfaction 
with  the  performance  of  existing  methods,  especially  originating 
from  inadequacies  that  arise  when  identification  of  the  microstruc¬ 
ture  of  geometrically  complex  regions  (e.g.,  fiber  tract  crossings)  is 
attempted.  One  reason  behind  the  inadequacy  is  the  creation  of  DW- 
MRI  models  of  Section  III  by  direct  adoption  of  DW-NMR  models. 

The  following  important  considerations  arise  in  the  course  of  the 
adoption: 

1.  In  contrast  to  DW-NMR,  in  MRI  the  signal  does  not  originate 
from  a  single  point  at  echo  time  (see  Section  II).  Even  in 
MRI  without  diffusion  weighting,  all  of  the  points  in  FID  fill 
up  the  k- space  and  therefore  all  of  them  are  necessary  to  con¬ 
struct  the  image  or  the  volume. 

2.  Imaging  gradients,  which  are  naturally  absent  in  NMR,  are 
usually  ignored  in  DW-MRI  models.  Attempts  to  integrate 
imaging  gradient  effects  have  unveiled  significant  problems 
in  the  derivation  of  DTI  model  (Ozcan,  2010a).  Moreover, 
these  problems  cannot  be  remedied  neither  with  optimization 
(Ozcan,  2011)  nor  with  nonlinear  estimation  (Ozcan,  2010b). 

3.  With  the  exception  of  Torrey’s  approach  (Torrey,  1956),  D  is 
treated  as  a  position  independent  quantity  within  all  of  the 
models.  For  example,  DTI  relies  on  the  assumption  that  “the 
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diffusion  coefficients  are  independent  of  space  (position) 
within  a  voxel”  [p.  253,  (Basser  et  al.,  1994)].  Otherwise, 
DTI  estimation  equations  cannot  be  derived.  Position  de¬ 
pendence  complicates  severely  the  solution  of  PDEs  for  DW 
imaging  models. 

4.  The  narrow  pulse  approximation  is  unrealistic.  Considering 
that  Wiener  process  has  nondifferentiable  continuous  sample 
paths,  the  irregularities  of  the  spin  motion  are  at  a  level  of 
complexity  which  is  unresolvable  at  any  time  scale  regardless 
how  short  are  the  duration  of  the  pulses. 

5.  With  the  exception  of  GDTI  proposed  by  Liu  et  al.  (2004), 
the  models  obtain  symmetric  quantities  either  by  projecting 
the  data  onto  symmetric  constructs  (e.g.,  ellipsoids  of  DTI), 
pruning  the  asymmetric  structures  (see  Section  III.B.  with 
spherical  harmonics  of  HARDI  or  odd  ranked  tensors  of 
DOT)  or  by  using  the  magnitude  of  the  signal  in  spectral 
methods  which  forces  the  Fourier  transform  to  result  in  (Her- 
mitian)  symmetric  outcomes.  An  appropriate  model  must  be 
unconstrained  so  that  any  property  such  as  symmetry  should 
be  revealed  by  the  data  themselves  establishing  an  evidence 
based  methodology.  Modeling  presumptions  constrain  the 
outcome  of  the  investigations. 

6.  Models  rely  on  Markovian  property  which  is  not  justifiable  in  re¬ 
stricted  environments  such  as  biological  tissue  (see  Section  II.A.). 

Motivated  by  these  concerns,  a  new  accurate  and  unifying  Fou¬ 
rier-based  theory  of  DW-MRI  has  been  developed  by  returning  to 
the  first  principles  of  DW-MR  signal  formation  in  (Ozcan,  2010c; 
Ozcan,  2011;  Ozcan  et  al.,  2011).  The  DW-MRI  signal  is  modeled 
by  expressing  the  evolution  of  the  phase  (Q/)  of  the  ith  spin’s  trans¬ 
verse  magnetization, 

rm(t)  =  exp  ( —jyClj )  m;  (lo )  (30) 

Eq.  (30)  is  the  Bloch  equation  obtained  when  relaxations  are 
neglected.  The  magnetization  changes  according  to  the  displace¬ 
ment  Wj  from  the  initial  position  xt  (t0): 

Xi(t)  =Xj(to)  +Wi(t),  (31) 

because  of  the  diffusion  sensitizing  magnetic  field  gradients  (see 
Section  I).  When  the  signal  formation  is  derived  for  the  PGSE 
experiment  (Ozcan,  2010c)  using  Eq.  (32)  and  (33),  the  initial  posi¬ 
tion  (t0)  and  the  ith  spin’s  displacement  integral,  Wf, 

rtdi  rtd4 

Wf  =  wMdx  -  /  Wi(x)dx,  (32) 

J  tdl  J  td3 

(see  Fig.  1  for  the  definition  of  tdi)  appear  in  the  formulations.  Fur¬ 
thermore,  the  derivations  in  (Ozcan,  2010c)  show  that  the  complex 
valued  DW-MRI  signal  that  comes  out  of  the  MR  scanner  is  the 
Fourier  transform  of  the  joint  distribution  of  the  number  of  spins 
starting  from  the  initial  position  x  and  possessing  W4  displacement 
integrals:  Pcfd  (x,  W4): 

Scfd(£mr  ,kD)=F[Pcfd(x,Wd)}  (33) 

Accordingly,  the  method  is  named  complete  Fourier  direct  (CFD) 
MRI  because  the  entire  (complex  valued)  data  constitute  the  Fourier 
transform  of  the  distribution  function  Fcfd  (x,  Wd).  In  Eq.  (33),  the 


Fourier  variables  kmv  and  kD  are  the  usual  MR  imaging  k-space  vari¬ 
able  (read-out  and  slice  select)  and  the  diffusion  sensitizing  gradient 
vector  &D  =  Gd  respectively.  Basically,  CFD-MRI  augments  the 
dimension  of  MR  imaging  by  adding,  on  top  of  position  information 
of  MRI  (see  also  the  discussion  on  paragraph  3  in  Section  I),  three 
more  dimensions  corresponding  to  displacement  integrals.  In  con¬ 
trast  to  the  existing  methods,  the  unifying  CFD-MR  framework 
does  not  separate  position  and  displacement  portions.  In  fact,  in  the 
derivations  of  CFD  model,  Pcfd  (x,  W4)  comes  up  naturally  as  a  joint 
distribution  function. 

It  is  very  important  to  note  that  Eq.  (33)  is  the  final  product  of 
intermediate  steps  that  guarantee  the  preservation  of  Scfd’s  Hermi- 
tian  symmetry  (Ozcan,  2010c).  This  property  is  the  only  physical 
evidence  at  hand  because  the  only  known  fact  is  that  Pcfd  is  a  real 
valued  function  being  the  count  of  (large  but)  finite  number  of 
spins.  This  implies  that  its  Fourier  transform  in  Eq.  (33)  must  be 
Hermitian, 

^cfd(^mr  5^d)  =  5cfd  (~kmx  ,  _ &d)  (34) 

(*  denotes  the  complex  conjugation). 

Rcfd  is  reconstructed  by  taking  the  inverse  Fourier  transform  of 
the  entire  data: 

Pea  (x,  Wd)  =  F-1  [Scfd  ( kml  ,  *D)] .  (35) 

CFD-MRI  addresses  the  concerns  raised  in  the  beginning  of  this 
section  as  follows: 

1.  CFD-MRI  calculations  incorporate  all  of  the  data.  The  signal 
is  not  calculated  from  a  single  point  but  from  all  the  points  of 
augmented  CFD  k-space  defined  by  &CFD  =  (k^,  &D).  £Cfd  is 
either  5  or  6  dimensional  depending  on  whether  slice  or  vol¬ 
ume  imaging  is  conducted.  The  last  3  dimensions  come  from 
3  dimensional  displacement  integral  vectors  defined  in  Eq. 
(32). 

2.  Imaging  gradients  are  included  by  means  of  kmr. 

3.  CFD-MRI  expresses  position  dependence  inherently  in  the 
first  argument  of  joint  distribution  function  Pcfd  (x,  W6)  which 
is  the  initial  position. 

4.  There  is  no  need  for  the  narrow  pulse  assumption.  In  Eq. 
(32),  the  integrals  are  evaluated  during  the  diffusion  sensitiz¬ 
ing  gradient  times. 

5.  The  signal  is  treated  as  a  complex  valued  quantity  in  contrast 
to  using  the  magnitude.  By  preserving  the  Hermitian  symme¬ 
try  during  the  evaluation  of  Eq.  (35),  asymmetry  is  allowed 
in  RCfd  (x,  W4).  This  guarantees  an  evidence  based  approach: 
if  there  is  asymmetry  or  symmetry  in  the  motion  characteris¬ 
tics,  data  will  show  it. 

6.  The  only  assumption  in  Eq.  (31)  is  the  continuity  of  the  dis¬ 
placement  Wi  since  a  spin  cannot  disappear  at  one  point  and 
reappear  at  a  different  one.  In  addition,  the  displacement  inte¬ 
gral,  W4,  is  not  a  stochastic  processes.  Therefore  any  property 
thereof,  such  as  being  Markovian,  is  not  needed  to  be  in  the 
model.  Accordingly,  CFD-MRI  points  out  that  the  propagator 
is  not  the  quantity  that  is  being  measured. 

The  simplicity  of  Eq.  (30)  (Bloch  equation)  and  Eq.  (31)  (the 
most  general  equation  of  motion)  make  their  combination  tractable 
by  use  of  particle  methods.  CFD-MRI  model  would  be  hard  to 


50  Vol.  22,  44-52  (2012) 


37 


Figure  3.  The  summary  of  the  DW  methods  is  shown  using  the  Fourier  transform  of  Fourier’s  portrait  on  top  right  corner.  MR  scanner  takes  the 
full  Fourier  transform  of  the  portrait.  CFD-MRI  reconstructs  the  portrait  by  using  the  inverse  transform.  Spectral  methods  take  the  magnitude  of 
the  transform  and  compute  the  inverse  transform.  The  outcome,  shown  in  the  middle  column  on  top,  does  not  resemble  the  original  portrait. 
Model  matching  methods  (arrow  at  the  bottom)  evaluate  the  inverse  transform  in  imaging  domain  and  work  in  a  mixed  physical-frequency 
space.  [Color  figure  can  be  viewed  in  the  online  issue,  which  is  available  at  wileyonlinelibrary.com.] 


derive  using  a  PDE  framework.  However,  theoretical  equivalence 
of  the  two  methods  (Feller,  1967)  indicates  that  particle  method 
derivations  are  sufficient  for  modeling. 

The  Fourier  relationship  established  with  CFD-MRI  clarifies  the 
mechanisms  of  existing  methods.  For  example,  the  Wiener  process 
that  describes  self-diffusion  in  a  liquid  has  displacements,  wh  that 
are  Gaussian  distributed  with  0  mean  and  covariance  matrix  equal 
to  the  diffusion  tensor  D.  It  can  be  derived  after  tedious  but  routine 
calculations  that  the  displacement  integrals,  Wf,  have  mean  equal 
to  0  and  covariance  equal  to  (Ozcan,  2011) 

E^Wf(Wf)Tj=?,2^A~^D  =  b,D.  (36) 

This  equation  shows  that  the  MR  scanner  acts  as  a  filter  on  the  dif¬ 
fusion  process  by  scaling  its  covariance  with  bt  that  only  depends 
the  duration  and  separation  times  of  the  diffusion  sensitizing  gra¬ 
dients.  The  DW  signal  that  comes  out  of  the  scanner  is  the  Fourier 
transform  of  the  distribution. 

Gaussians  are  eigenfunctions  of  the  Fourier  transform,  i.e.  the 
Fourier  transform  of  a  Gaussian  is  also  a  Gaussian  (Ozcan,  2011): 

F[exp ((W?)T(b,D)-lW?)\  ssexp {{kD)T{b,D)kD).  (37) 

Equation  (37)  is  key  in  understanding  the  model  matching  methods. 
Equation  (9)  is  obtained  by  evaluating  the  Fourier  transform  in  half  of 
the  variables,  namely  the  imaging  portion  defined  by  kmr  and  leaving 
the  displacement  integral  portion  kD  intact  In  this  mixed  physical-fre¬ 
quency  variable  setup,  the  displacement  portion  remains  on  the  Fourier 
domain.  According  to  the  right  hand  side  of  Eq.  (37),  DTI  estimates  D 


in  Fourier  domain  more  easily  than  in  the  physical  domain  on  the  left 
side  which  involves  of  the  inverse  of  D.  In  more  complex  microstruc¬ 
tures,  high  order  approximations  follow  the  same  strategy.  They  attempt 
to  expand  the  mathematical  Fourier  transform  by  adding  higher  order 
terms  described  in  Section  9  to  the  argument  of  the  eigenfunction. 

According  to  Eq.  (35),  CFD-MRI  reconstructs  Pcfd  using  the 
inverse  Fourier  transform  (discrete  Fourier  transform  is  used  in 
practice).  Spectral  methods  of  Section  IV  are  based  on  the  same 
principle  but  CFD  evaluates  the  Fourier  transform  correctly  (see 
item  5  above),  i.e.,  using  complex  valued  signals.  It  obtains  the  cor¬ 
rect  distribution  function  (see  item  6  above),  in  the  sense  that  Fcfd 
being  a  joint  distribution,  cannot  necessarily  be  factored  into  a  spin 
density  and  displacement  descriptor.  The  mechanisms  behind  the 
methods  are  summarized  in  Figure  3. 

V.  CONCLUSIONS 

All  DW-MRI  models  process  the  data  with  the  same  goal  in  mind: 
noninvasive  description  of  micro  structural  changes  that  will  endow 
extraction  of  relevant  information  pertaining  to  diagnosis,  prognosis 
and  management  of  important  health  problems  from  different  areas. 
The  large  number  of  models  described  herein  is  proof  of  continuing 
research  effort  in  the  quest  of  a  more  adequate  model  for  DW-MRI. 
The  unification  of  the  plethora  of  approaches  will  improve  the  accu¬ 
racy  of  the  assessment  provided  by  the  modality.  In  this  manuscript, 
the  methods  are  summarized  by  describing  their  mathematical  prop¬ 
erties.  In  addition,  they  are  interpreted  using  the  Fourier  based 
CFD-MRI  method  with  the  aim  of  reaching  a  consensus  among  the 
methods.  An  important  goal  is  the  incorporation  of  these  new  find¬ 
ings  into  the  tractography  methods,  which  is  reserved  for  future 
studies. 
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Appendix  4 

Importance  of  Sleep  to  Humans 
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WHAT  HAPPENS  TO  THE  BRAIN  WHEN  WE  SLEEP? 


Brain  electrical  activity  changes  in  a  highly  stereotypic  pattern 


Wake 
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WHAT  HAPPENS  TO  THE  BRAIN  WHEN  WE  SLEEP? 


The  brain  reduces  its  energy  use 
in  deep  sleep 
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Maquet  (1995),  Braun  et  al.  (1997) 


This  reduced  energy 
consumption  is  greatest  in  parts 
of  the  brain  involved  in  attention, 
cognition  and  integration  of 
information 


Frontal 


Parietal 


Nofzinger  (2002) 


Temporal 


The  transfer  of  information 
between  brain  regions  -  the 
functional  connections  between 
brain  networks  -  becomes  less 

global  and  more  local  (Horovitzet 
al.,  2009;  Spoormaker  et  al.,  2010;  Larson- 
Prior  et  al.,  2009,  2011) 
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WHY  DO  WE  SLEEP? 


We  do  not  know  WHY  we  sleep,  but  we  know  that  it  is  important,  because 


All  organisms  sleep: 


Even  if  it  is  life  threatening! 
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WHY  DO  WE  SLEEP? 

And  when  we  don’t  sleep: 

Hours  sleep/night  Daytime  effect 


Possible  long-term 
consequences 
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Performance  improvement 


WHAT  ARE  THE  BENEFITS  OF  SLEEP? 


Sleep  promotes  at  least  some  types  of  learning  and  memory 


Sleep  improves  memory  of  a  learned  sequence  of 
finger  taps  to  a  greater  extent  than  is  seen  with  the 
passage  of  time  alone 


Walker  and  Stickgold  (2006)  Ann  Rev  Psychology  57:139 


No  sleep,  day  1 


After  sleep,  day  2 
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WHAT  ARE  THE  BENEFITS  OF  SLEEP? 

Sleep  promotes  good  mental  and  physical  health 


HYPOSOMNIA 

I 

Sleep  Fragmentation 

I 


Cognitive  Dysfunction 
Executive  Dysfu 


HYPERSOMNIA 

Increased  risk  of  obesity 

Metabolic  Syndrome 

4 


MAJOR  DEPRESSION 


STROKE 
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HOW  MUCH  SLEEP  DO  WE  NEED? 


Humans  sleep  between  4  and  10  hours/night  -  there  is  no  number  of  hours 
that  is  “right”  for  everyone 


The  amount  of  time  we  need  to  sleep  changes  with  age 


How  Much  Sleep  Do  You  Really  Need? 

Age 

Sleep  Needs 

Newborns  (0-2  months) 

12-18  hours 

Infants  (3  to  11  months) 

14  to  15  hours 

Toddlers  (1-3  years) 

12  to  14  hours 

Preschoolers  (3-5  years) 

11  to  13  hours 

School-age  children  (5-10  years) 

10  to  11  hours 

Teens  (10-17) 

8.5-9.25  hours 

Adults 

7-9  hours 

Source:  National  Sleep  Foundation 
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